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' No surface is perfectly planar at all scales. The notion of flatness of a surface therefore depends on 

I the size of the probe used to observe it. As a consequence rough interfaces are abundant in nature. 

, Here the old, but still active field of rough surface scattering of electromagnetic waves is addressed. 

2f^' This topic has implications and practical applications in fields as diverse as observational astronomy 

^ ' and the electronics industry. This article reviews the theoretical and computational foundation and 

, methods used in the study of rough surface scattering. Furthermore, it presents and explains the 

I ■ physical origin of a series of multiple scattering surface phenomena. In particular what is discussed 
are: the enhanced backscattering and satellite peak phenomena, coherent effects in angular intensity 

I— I, correlation functions and second harmonic generated light (a non-linear effect). 
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I. INTRODUCTION 

We are surrounded by waves, and they effect our daily life in a way that many of us are not aware. Sound and light 
are our main tools for observing our immediate surroundings. Light is for example responsible for you being able to 
read these lines, and more important, to get access to the vast majority of the knowledge accumulated in writings 
by man throughout centuries of intellectual activities. X-ray and ultra sound techniques have given tremendous 
contribution to the success of modern medicine. Radio- and micro-waves are invaluable in modern communication 
technology including cellular phones and radio and TV broadcastings. Understanding of quantum waves, and their 
behavior, constitutes the foundation of electronics and semiconductor technologies — an essential ingredient in the 
past and future progress of computer hardware. The above list is not at all, or intended to be, complete. It could in 
fact easily been made much longer. However, the bottom line that we want to make here is that with the ubiquitous 
presence of wave phenomena in various applications, it is not surprising to find that wave phenomena have had, and 
still have, a prominent position in our studies of the physical world, and even today such phenomena are of out-most 
importance in science, medicine and technology. 
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FIG. 1: An illustration showing the transition from specular (Fig. la) to diffuse scattering (Fig. Id) of light from a surface of 
increasing roughness. The arrows indicate the direction of the incident light (After Ref. [5]). 



If you take an average introductionary text on wave phenomena, you will find discussions of how plane waves 
of constant frequency propagates in a homogeneous, isotropic medium. Thereafter, the authors typically discuss the 
scattering and transmission of such waves at a planar interface separating two semi-infinite media of different dielectric 
properties [11] — the Fresnel formulae. These formulae serve to accurately describe the scattering of light from for 
example a mirror. However, from our everyday experience, we know that most surfaces are not mirror like, and 
naturally occurring objects are more complicated then two semi-infinite media. Most naturally occurring surfaces are 
actually not smooth at all. They are, however, rough in some sense. In fact, all objects, man-made or not, have to be 
rough at atomic scales, but such small length scales are normally not resolved by our probes. 

It should be kept in mind that the characterization of a surface as rough or smooth is noticeably not unique, and 
it is not a intrinsic property of the surface. Instead, however, it depends on the wavelength used to "observe" the 
surface. If the typical roughness is on a scale much smaller then the wavelength of the probe, this surface is considered 
as smooth. However, by reducing the wavelength of the light, the same surface might also be characterized as being 
rough. It is, among other factors, the surface topography and the wavelength of the probe, as we will see below, that 
together go into the characterization of a surface as being rougli^ . 

Let us from now on assume an electromagnetic probe, i.e. light. If the surface can be considered as smooth, 
light is scattered (coherently) into the specular direction. As the roughness of the surface is increased so that the 
surface becomes weakly rough, a small fraction of the incident light will be scattered into other directions than the 
specular one. This non-specular scattering is called diffuse scattering or by some authors incoherent scattering. As 
the roughness is increased even further, the diffuse (incoherent) component of the scattered light is increased on the 
expense of the specular component. When the surface roughness is so that the specular component can be more-or-less 
neglected as compared to the diffuse component, the surface is said to be strongly rough. This transmission from as 
smooth to a strongly rough surface is depict in Figs. 1. 

Due to the practical applications of waves, and the number of naturally occurring surfaces being rough, it is rather 
remarkable that it took several hundreds years from the birth of optics as a scientific discipline to someone started to 
considered wave scattering from rough surfaces. As far as we know today, the first such theoretical study was made 
at the end of the 19th century (probably in the year of 1877) by one of the greatest scientists of its time, the British 
physicist Lord Rayleigh [1, 2]. He considered the scattering of light incident normally onto a sinusoidal surface. 

In 1913 Mandel'shtam studied how light was scattered from liquid surfaces [3]. By doing so, he became the first 
to consider scattering from randomly rough surfaces. This, as it turned out. should define the beginning of an active 
research area wave scattering from randomly rough surfaces — which still today is an active field. However, it was 
first after the last world war that the research effort put into the field stated to accelerate [4]. Since that time, a 
massive body of research literature has been generated in the field [5-8] . 

Up to the mid 1980's most of the theories used in this field were single scattering theories [5-7]. However, from 
then on the main focus of the research has been on multiple scattering theories. In addition, advances in experi- 
mental techniques has lately enabled experimentalists to fabricate surfaces under well controlled conditions by using 
a holographic grating technique [9]. This has opened up a unique possibility for direct comparison of theory and 
experiments in a way not possible a few decades ago. 

Inspired by the works of Lord Rayleigh [1, 2] researchers developed a criterion — the Rayleigh criterion — that 



^ When discussing the Rayleigh criterion later in this section we will see that also the angle of incidence of the light will play an important 
role. 
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could be used to determine when a given surface was to be considered as rough. Here both the wavelength of the 
incident light as well as its angle of incidence are incorporated [6] . 

To illustrate how this comes about, let us consider a rough surface defined by xj, = ((xi). On this surface we pick 
two arbitrary points C(0) ^'^d (C) C(C'))- could now be asked: What is the phase difference between two waves 
being scattered from these two points? For simplicity we will here only consider the specular direction. Under this 
assumption it is straight forward to show that the phase difference is given by the following expression 



where |k| — 2tt/X is the modulus of the wave vector of the incident light of wavelength A, and is the angle of 
incidence of the light as measured from the normal to the mean surface. From Eq. (1) we immediately observe that if 
the surface is planar, so that = CiO^ the phase difference (in the specular direction) is always zero independent 
of the angle of incidence. However, if the surface is rough, 7^ in general. If A(f) <C tt, the two waves will be in, or 
almost in, phase and they will thus interfere constructively. On the other hand, if A<f> ~ tt, they will be (morc-or-lcss) 
completely out off phase and as a result interfere destructively, and no, or almost no, energy will be scattered into 
the specular direction. In terms of the phase, a smooth surface would correspond to <^ tt, and a rough one to 
A(f) ~ TT. Thus, A(p = 7r/2 might be considered as the borderline between a smooth and a rough surface; if A(f) < 7r/2 
the surface is smooth, and otherwise (7r/2 < A(j) < tt) it is rough. The criterion A(j) < 7r/2 is the famous Rayleigh 
criterion for a smooth surface. 

If the surface is randomly rough, it is practical to replace the height difference — C(C') by a typical height 
fluctuation as provided, for example, by the RMS-height, 6, of the surface. Hence, the Rayleigh criterion can be 
expressed as 



where Ra is the so-called Rayleigh parameter. From the Rayleigh criterion, Ra < 7r/4, it should be observed that 
in addition to the surface topography itself and the wavelength of the light, also its angle of incidence goes into 
determining if a surface is rough or not. This is probably the most important lesson to be learned today from the 
Rayleigh criterion. 

The present review consists of basically two main parts — one focus theoretical methods whilst the other one is 
devoted to rough surface scattering phenomenology. In the first part we try to present an overview of some of the 
main theories and methods used in the study of wave scattering from randomly rough surfaces. We start in Sect. II by 
recapitulating the basic results of electromagnetic theory including Maxwell's equations. This section serves among 
other things to define our notation. Then we continue by describing how to characterize randomly rough surfaces 
(Sect. HE). Sect. HI is devoted to the quantities and main techniques used in the field of electromagnetic wave 
scattering from randomly rough surfaces. We here review classical theories like small amplitude perturbation theory, 
many-body perturbation theory as well as numerical simulation approaches. Finally in Sect. IV we discuss some of 
the phenomena that may occur when light is scattered from rough surfaces. Such effects include the backscattering 
and satellite peaks phenomena (weak localization), Anderson localization, angular intensity correlation effects and 
nonlinear effects (second harmonic generation). 



The present review mainly concern itself with rough surfaces and the scattering of electromagnetic wave from such. 
In this section we therefor review some of the basic results of electromagnetic theory, including surface polaritons. 
The present section also serves to define our notation that we will use extensively in the following sections. The style 
of this review is kept quite brief, since all the material should be well known. A more thorough treatments can be 
found for example in the classical text on electrodynamics by J. D. Jackson [34]. 



A<p = 2|k||C(0-C(r)|cos^o, 



(1) 




(2) 



II. ELEMENTS OF ELECTROMAGNETIC THEORY 



A. 



Maxwell'i 



's Equations and the Constitutive Equations 



1. Maxwell's Equations 



The Maxwell's equations, which unify in one magnificent theory all the phenomena of electricity and magnetism, 
were put forward by the Scottish physicist James Clerk Maxwell (1831-1879). These equations are the fundamental 



5 



Quantity 


Sl-unit 


Name 


E 


V/m 


Electric field 


H 


A/m 


Magnetic field 


D 


C/m'' 


Electric displacement 


B 


Wh/m^ 


Magnetic induction 


P 


C/m^ 


Charge density 


J 




Current density 



TABLE I: Summary of the quantities contained in Maxwell's equations, as well as their Sl-units. 



equations of electromagnetism, in the same way that Newton's law is to classical mechanics. In fact, the Maxwell's 
equations are in a way even more fundamental since they are consistent with the theory of special relativity that 
Einstein develop years later. Because all of electromagnetism is contained within this set of equations, they are 
definitely among one of the greatest triumphs of the human mind. 

Strictly speaking the equations put forward by Maxweh only applies to point charges in vacuum. A dielectric, for 
example, is a collection of a very huge number of point charges. To deal with them all individually is an impossible 
task. It is therefore practical to introduce effective fields, D and H, to represent their collective behavior. This 
dielectric approach to electromagnetism represents great simplifications for many (near-to-natural) systems. It is 
based on the following two assumptions [34]: (i) the response of the background medium is dipole like as well as 
linear in the applied fields, and {n) the medium is homogeneous (or close to) throughout a given region. The first 
assumption obviously breaks down if the fields becomes to strong while the latter breaks down on short length scales. 
Hence the resulting effective field theory, or effective Maxwell theory as we might call it, should be treated as a long- 
wavelength approximation to electromagnetism for weak fields. In most practical situations the above approximations 
are fortunately well satisfied, and in particular they are valid for the type of scattering system that we will be 
considering. 

In the Sl-system, the (effective) Maxwell's equations take on the following form: 

V-D = p, (3a) 
V-B = 0, (3b) 

V X E = - — , (3c) 

r\T~V 

V X H = + J- (3d) 

Here E and H denote the electric and magnetic field vectors respectively. These field vectors make together up what 
is known as the electromagnetic field. The field quantities D and B, known as the electrical displacement and the 
magnetic induction respectively, are included in order to describe the effect of the electromagnetic field on matter. 
Finally, p and J denote the charge density and the current density respectively. Those two latter quantities act like 
sources for the electromagnetic field, E and H, and they fulfill the continuity relation 

|+V.J = 0. (4) 

The various quantities appearing in the Maxwell's equations, and related formulae, are summarized in Table HA 
where also their Sl-units are given. 



2. Constitutive Equations 

The Maxwell's equations (3) consist of eight scalar equations. However, on the other hand the field vectors, E, H, 
D, and B, represent in total 12 (scalar) variables, 3 for each of the 4 vectors. Thus, obviously, the Maxwell's equations 
alone do not uniquely specify a solution. Therefore, in order to obtain a unique solution to the Maxwell's equations, 
those are supplemented by so-called constitutive relations also known as material equations. These relations read 



D = eE, 
B = /xH. 



(5a) 
(5b) 
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FIG. 2: A sketch of a general interface separating two dielectric media. 



Here e and ji are the constitutive parameters which are tensors of 2nd order and known as the permittivity^ , and the 
permeM,hility tensor respectively. In general these tensors are rather complicated functions of the spatial variable x 
and the field vectors E and H. However, for an isotropic and homogeneous medium, these tensors reduce to scalars, 
and if the field-strengths are not too large, they may be considered as independent of the field vectors. In this latter 
case we arc dealing with linear electromagnetic theory. The fascinating, but complicated nonlinear electromagnetic 
theory [35] where per definition e and /x depend on E and H, will not be discussed here in any depth. 
Eqs. (5) can within linear electromagnetic theory be cast into the equivalent form 

D = £oE + P, (6a) 
B = /ioH + M, (6b) 

where P and M are the electric and magnetic polarizations respectively. The constants so and /Uq are the permittivity 
and permeability of vacuum respectively. In the Sl-system they have the following values 

£0 = 8.854 X 10-^2 -p/ra, (7a) 
Mo = 47r X 10"^ H/m. (7b) 



B. The Electromagnetic Wave Equations 

Probably the two most important consequences of the Maxwell's equations are the wave equations and the existence 
of solutions to these which are known as electrom,agnetic waves due to the wave-like nature of such solutions. In this 
section we derive the wave equations in a material medium. For simplicity, and since it is the most relevant case for 
this review, we will limit ourselves to a region of space which is source free and isotropic. 

The derivation of the wave equations for the E-field in a source- free region (i.e. p = and J = 0), is achieved 
by ehminating the H-field from the Maxwell's equations. This is done by taking the curl of Eq. (3c), substituting 
Eq. (3rf), and taking advantage of the constitutive relations (5). The result is 

9^E 

V X (V X E) + £/x-^ = 0. (8) 

By applying the vector identity V x (V x A) = V(V • A) — V^A to Eq. (8) and taking advantage of Eq. (3a) we 
obtain the well-known standard (space-time) wave equation for the electrical field in a source-free, homogeneous and 
isotropic medium 

V2E-/.£^ = 0. (9) 

In a similar way one can obtain a wave equation for the magnetic field by eliminating the electric field from the 
Maxwell's equations. 



^ This quantity is also known as tiie dielectric function. 
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It should be notice that not every solution to the wave equation is also a solution to the Maxwell's equations. For it 
to be, it must in addition satisfy Gauss's law, V-E = k- E = Oin order to also be a solution of Maxwell's equations^. 
As the reader readily may check the wave equation has a solution E = exp(ik • r — iuit) if a; = ck. This solution is the 
plane wave solution. 



C. Boundary Conditions 

In Sect. II A we introduced the Maxwell's equations and the constitutive relations. These equations can be solved 
for the field vectors in a region of space containing no boundaries. However, no real media are infinite, i.e. without 
boundaries. For practical applications of the electromagnetic theory it is therefore important to know how to treat 
the boundaries between two media of different electromagnetic properties. It is this question that we will address in 
this section. 

Let us consider the geometry of Figure 2. It shows an arbitrary interface separating the otherwise homogeneous, 

isotropic and linear media labeled ±. We have also introduced a normal vector for the interface, n, which is directed 
into medium +. The question we now address is: How are the electromagnetic field vectors for the two media in 
the immediate vicinity of the interface related to each other? The answer to this question should be well-known and 
can be found in nearly any book on electromagnetic theory, e.g. in Refs. [11] and [34]. The results, for which the 
derivation will not be repeated here, are 



n (B_ 


-B+) 


= 0, 


(10a) 


n (D_ 


-D+) 


= Pa, 


(10b) 


nx (E_ 


-E+) 


= 0, 


(10c) 


nx (H_ 


-H+) 




(lOd) 



where the vector subscripts, ±, are referring to the media where the field vectors are evaluated. In Eqs. (10) and 
3s denote the surface charge density and the surface current density respectively, while the other quantities have been 
defined earlier. In many areas of optics one deals with situations where the surface charge density and the surface 
current density are zero. Under such circumstances the normal component of B and D are continuous, while the 
vectors E and H have continuous tangential components. 

It should be stressed that in arriving at the results (10), it has been assumed that the electromagnetic properties 
take on their bulk values all the way to the surface. This is obviously not true, but is a good approximation whenever 
the mean field theory applies. 



1. Boundary Condition at a General One-Dimensional Surface 



Most of this review will concern itself with randomly rough surfaces that are effectively one-dimensional, i.e. the 
surface profile function C, has a non-trivial dependence only on xi, say, and does not depend explicitly on X2. In this 
case the boundary conditions (10) simplifies somewhat. This is what we plan to outline in this section. 

Let us start by assuming, without loss of generality, that the plane of incidence is the xix^-plane and that the 
incident light is either p- or s-polarized. In such case, there is only one non-trivial field component needed in order 
to fully describe the electromagnetic filed. For p-polarization this component is H2, while for s-polarization it is E2. 
Thus the primary field for a one- dimensional interface problem can be written as 

$.(.„.3M = (ff^'^f)' ^ = ^' (11) 
[ i;2(a;i,a;3|a;), u = s, 

where a harmonic time-dependence, exp(— iwi), has been assumed, but suppressed. This form for the primary field 
will be used frequently throughout this review. Notice the fact that the primary field can be fully described by a 
single vector component. This represents a dramatic simplification of the problem since it is reduced form a vector 
problem down to a scalar one. 



^ One explicate example of this is provided by E = 63 Eq cos(fea;3 — tot) that satisfies the wave-equation, but not V • E = 0. It must 
therefore be discarded as a solution of the Maxwell's equations. 
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When <I>i^(.xi, x^luj) is known, the remaining component of the electromagnetic field can be calculated from it alone. 
These components are given for p-polarization by 

i d 

Ei{xi,X3\u>) = — H2{xi,X3\u>), (12a) 

i d 

^3(a;i,a;3|a;) = — --^-—H2{xi,X3\u}), (12b) 
uje(uj)oxi 

and for s-polarization 

t d 

Hi{xi,X3\lo) = -— - — E2{xi,X3\uj), (13a) 

i 

H3{xi,X3\uj) = -— E2{xi,X3\(j). (13b) 

uj^[ijj) 0x1 

In the above equations s{lo) and fi{Lu) denoted the dielectric function and the magnetic permeability respectively of 
the medium where the fields are being evaluated. The relations (12) and (13) arc easily derived by using the two 
curl-equations contained in the Maxwell's equations, Eqs. (3c) and (3d), together with the constitutive relations, 
Eqs. (6). 

Let us now try to focus on the boundary conditions that the primary field 4>^(xi, xsjo;) will be subjected to. By 
construction ^^{xi,X3\lj) is a tangential field independent of polarization. Therefore it follows automatically from 
Eqs. (10a) and (lOd) {ps = = 0) that 

^t{xi,X3\w)\^^^^^^^^ = $-(xi,X3|u;)|^^^^(^^^, (14a) 

where C,{^i) denotes the interface separating the two materials of different dielectric properties. 

In order to satisfy the remaining boundary conditions expressed in Eqs. (10), we notice that for respectively p- and 
s-polarization we have 

1 

n X E = 62 -j—dn^p, 

n X H = 62 — ^-—dn^s, 

where 9„ denotes the normal derivative to the surface. If the one-dimensional interface can be represented as X3 = 
C,{xi), where Ci^i) is a single- valued function of x\ the normal derivative becomes 

Vi + (C'(2^i))^ 

where dx^ = d/dxi and 

^ C'(a;i)ei +63 

Here dj are the standard unit vectors. Hence the remaining boundary conditions can be expressed as 



(14b) 



^ dr,^t{XuX3\u^)\.,. = — dn^. (^1 ' ^3 1 w) I . , (14c) 



where kJ(w) are defined as 



«J(-) = r^^"^; (i4d) 

[ At±(a;), v = s 

Eqs. (14) are the final result for the boundary conditions to be satisfied by the primary field (^i,{x\,X3\uj) on a 
one-dimensional interface X3 = C,{xx). 
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D. Surface Plasmon Polaritons 



In subsequent sections, we will see that so-called surface plasmon polaritons, or for short just SPPs, will play an 
important roll for the rough surface scattering problem. We will therefore in this section define and discuss some of 
the distinguishing properties of such modes. 

Before starting our discussion, we have to know what a polariton is: According to its classical definition a polariton 

is defined to be an elementary electromagnetic wave, and therefore a solution of the Maxwell's equations, that may 
couple to one of several possible excitations possible in a condensed medium. Examples of such excitations are 
plasmons, phonons, mognons etc., and in such cases one talks of plasmon polaritons, phonon polaritons and magnon 
polaritons. With the notion of polariton established, one might definition an SPP as follows: A surface plasmon 
polariton is a plasmon polariton where the associated electromagnetic field is confined to the surface separating two 
dielectric medium. 



1. SPPs on a plan surface geometry 

To see under which condition SPPs might exist, and to discuss some of their properties, we will consider a planar 

interface separating two isotropic and homogeneous media. For simplicity, the coordinate system will be chosen so 
that the interface is located at = 0. The materials above {x^ > 0) and below {x^ < 0) this surface will be 
characterized by frequency dependent dielectric functions e+ (w) and S- (w) respectively. For simplistic reasons, which 
are not essential for the present discussion, we will assume that the imaginary part of the dielectric hmctions can be 
neglected. The conclusion that we arrive at herein will, however, be independent of this assumption. Furthermore, 
we will assume either pure p- or s-polarization of the incident light. Hence the scalar wave equation might be used. 
A more complete discussion using vector fields can be found in Refs. [36] and [37] . 

According to the definition of SPP, we are interested in solutions to the Maxwell equations, equivalent in our case 
to the scalar wave equation, that are wave-like parallel to the surface x^ = and that decays exponentially with 
increasing distance from the surface into each of the two media. Such a solution can be represented as 

^tixi,X3\to) = ^±eT&(-)-3g«fcxi^ u=p,s, (15) 
where represents the amplitudes (to be determined). The decay constants l3±{co) are defined as 



/?±(w) = (16) 



and they must be real and positive for Eq. (15) to describe an electromagnetic wave localized to the surface . To 
investigate if Eq. (15) is an acceptable solution for our scattering system, we have to impose the boundary conditions, 
given in Eqs. (14), for the two polarizations. By utilizing the continuity of the fields on the flat surface, Eq. (14a), 
one finds that 

for all locations along the surface. Moreover, the normal derivative condition, Eq. (14c), gives the following condition 
for the existence for surface plasmon polaritons on a flat surface (9„ = dx^ ) 

M-K!LM)a. = o. (17, 



where we recall the definitions of k^{uj) from Eq. (14d). The most immediate consequence of this relation is the 
following: Since for s-polarization Kf{uj) = /xq (for non-magnetic materials) and /3±(a;) are assumed to be real and 
positive, the only solution to Eq. (17) is As = 0. Thus, a s-polarized surface plasmon polariton (surface wave) cannot 
exist for the scattering system that we are considering. 

However, for p-polarization, where k^{oj) = e±(w) a non- trivial solutions might exist. It is given by (assuming that 
Api^Q) 

£-kM _ _^+M /.ox 



* If we had allowed the dielectric functions of the problem to be complex with Ime±{w) > 0, we would have to require that Rel3±{ui) > 0. 



10 



which is the dispersion relation for surface plasmon polaritons on a flat interface. For this relation to be satisfied, since 
f3±{(jj) are both assumed to be positive, the two dielectric functions of the scattering system have to have different 
signs due to the presence of the negative sign on the right-hand-side of Eq. (18). Only such combination of materials 
will support surface plasmon polaritons. An important example of such a system at optical frequencies, is a metal 
with a planar interface to vacuum. 

By squaring both sides of Eq. (18) as well as taking advantage of Eq. (16), the dispersion relation can be expressed 

as 



ksp{u>) = ±\ — --— r^ — - (19) 

This equation gives an explicit expression for the wave vector of the surface plasmon polariton. However, this formulae 
should be used with some care since it may, from the way it is derived from Eq. (18), introduce some spurious solutions. 
The additional and sufficient requirement that have to be satisfied is that f3± {w) are positive while the two dielectric 
functions, s±{w), have different sign. 



2. SPPs at a planar free electron metal surface 

Let us for illustrative purposes consider a free electron metal with a planar interface to vacuum. For such a metal 
the dielectric function is known to be [38] 



(20) 



while the one for vacuum is e+(w) = 1- In the above equation £00 (w) is the background dielectric constant of the 
material while Wp is the electronic plasma frequency. With Eq. (20) the frequency of the SPP can be shown to be 




It should be noticed from this equation that 



_1 



(21) 



kc, 



£00 (w) + 1 



k 
k 



0, 



00, 



(22) 



This means that in the small wave vector limit the surface plasmon polariton is photon-like, while it in the large 
wave vector limit it is plasmon-like. In Fig. 3 the dispersion relation, Eq. (21), for a free electron metal is plotted. In 
this figure we have also included the light-line (dash-dotted line) iv = kc, as well as the large wave vector limit (dashed 
line) of ujsp{k). 

From Fig. 3 we see that the dispersion curve for the SPP lies entirely to the right of the light-line, uj = kc. The 
physical consequence of this is that there is no coupling between the surface plasmon polariton and light in vacuum 
for a flat vacuum-metal interface. Or put another way, light incident onto a planar vacuum-metal interface cannot 
excite surface plasmon polaritons. Later, however, we will see that if the surface is rough such coupling is possible. 
This will give rise to many new and interesting multiple-scattering effects, as we discuss in some detail in Sect. IV. 



E. Characterization of Random Rough Surfaces 

Almost everyone grows up with some kind of intuitive "feeling" of what is meant by a rough surface. A fractured 

stone, say, is normally looked upon as being rough, while a piece of paper as being smooth. However, on the micro- 
scale, where the human eye is not very sensitive, also the paper has some kind of structure. So in a strict sense, 
both the paper and the stone surface are rough. Paper is made out of flbers which is quite different from the crystals 
that are seen on the micro scale of the surface of the fractured stone. So the question is: How shall we quantify the 
difference in roughness between say the paper and the stone surface? One possibility is to measure by some suitable 



11 



0.8 



/ 




0.6 



3"" 0.4 
3 



0.2 











0.5 



1 



1.5 



2 



2.5 



3 



p 



FIG. 3: Dispersion relation curve, Eq. (21), for a surface plasrnon polariton (solid line) at a flat interface between a simple 
metal and vacuum (assuming £oo(<x') = 1). The dash-dotted line represents the light line lo = kc, while the dashed line is the 
large momentum asymptotic limit cJsp = uJp/V2 (see Eq. (22)). 



technique the surface topography. Such measurements will of course produce different results for the paper and stone 
surface. However, if we move to another area of the fractured stone and measure the surface topography here, we 
will obviously not get the same result as obtained in the previous measurement taken from another area of the same 
surface. So, how shall we be able to characterize the rough surfaces at hand, so that we are able to distinguish them 
from each other? In this section we intend to discuss in some detail how to characterize randomly rough surfaces in 
a quantitative way. 

However, before we do so, let us take a look at what kind of rough surfaces we have. Depending on how the surface 
height fluctuates around some reference surface, we may categorize them as being deterministic or randomly rough. 
For random surfaces, one may in addition group them as correlated or uncorrelated surfaces, and they might occur 
as fractal or non-fractal surfaces depending on under which conditions they were formed. Rough surfaces that are 
found in nature are normally randomly rough, correlated surfaces. We will therefore proceed by discussing how to 
characterize such surfaces. 



Two randomly rough surfaces arc never identical. Thus the knowledge of the surface topography alone is therefore 
not enough to be able to say if two rough surfaces were generated by the same underlying process, and therefore have 
to be looked upon as being identical. However, if we assume that the randomly rough surface can be considered as a 
continuous random process [12-17], then a statistical description might be relevant and useful. We will now introduce 
this method of characterization. 

Under experimental conditions, the surface topography is measured relative to some reference surface. In our case 
we will assume that this reference surface is a planer surface. Other choices might be practical in some cases, but this 
will not be discussed here. Furthermore, it is convenient to choose our coordinate system so that this planar surface 
is located at X3 = 0. In this case the randomly rough surface is just the roughness that perturb the plane X3 = 0. For 
simplicity, we limit our discussion to one-dimensional surfaces. The extension to (isotropic) two-dimensional surfaces 
is trivial. For the purpose of this introduction it will be assumed that the surface does not possess any overhangs^, 
that is to say that the surface profile function, that we will denote by C(2;i), is a single- valued function of the lateral 



F. A Statistical description of Randomly Rough Surfaces 



Such surfaces are also known as reentrant surfaces. 
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coordinate xi. For characterization of surfaces where the surface profile function does not fulfill this property the 
reader is invited to consult Ref. [18]. 

In order for the surface profile function C(^ci) to be planar on average, it must, with our choice for the coordinate 
system, have a vanishing mean, i.e. we must require that 

(C(xi)) = 0. (23) 

Here the angle brackets are used to denote a spatial average over a large spatial region. If, however, the surface is 
ergodic [13-17], as we will assume here, this spatial average is equal to an average over an ensemble of realizations of 
^(xi). It is therefore, imder the assumption of ergodicity, more convenient to think of (•) as an ensemble average. 

Another notion that is important when characterizing rough surfaces, and studying light scattering from such, is 
the one of stationarity [6]. A surface is said to be stationary, or translation invariant, if its statistical properties are 
independent of which portion of the surface was used in their determination. That the surface roughness possess 
stationarity is crucial for the applicability of many of the theories used to study rough surface scattering. Rigorous 
numerical simulations (Sect. Ill I), however, can still handle non-stationary surfaces. 



1. Gaussian Random Surfaces 

In theoretical studies of light scattering from rough surfaces, the random surfaces have in the overall majority of 
the studies been assumed to possess Gaussian height statistics. Such a statistics is rather appealing from a theoretical 
point of view since all moments can be related to the two first moments, such moments either vanish (odd moments), 
or they are related to the second moment (even moments) [6] . 

The zero- mean property, Eq. (23), does not specify how the different heights are located relative to one another 
along the surface. Such information is provided by the height-height correlation function. Under the assumption of 
C{xi) being stationary we can write 

(C(a;i)C(x;)) = 5'W{\x,-x[\), (24) 

where 5 is the RMS- height of the surface profile function, and VF(|xi|) is the height auto-correlation function normalized 
so that W{0) = 1. In cases where iy(|a;i|) = 1 (H^(|a;i|) = —1) one speaks of perfect correlation (anti-correlation). 
Furthermore it can be shown that — 1 < W^(|xi |) < 1. Notice, that since the heights-distribution is Gaussian, Eqs. (23) 
and (24) together determines uniquely the statistical properties of the surface since all higher order moments can (for 
a Gaussian surface) be related to the first two. 

In many of the perturbation theories developed for rough surface scattering, the power spectrum of the surface 
randomness is a quantity that appear more-or-less naturally. It is defined as the Fourier transform of the (normalized) 
correlation function 

/oo 
dx,W{\x,\)e-''''\ (25) 
-oo 

In order to get an intuitive picture of how the surface height varies along the surface, it is often useful to supply 
the mean slope, s, and the mean distance between consecutive peaks and valleys, {D), as measured along the (lateral) 
a;i-direction. For a stationary zero-mean, Gaussian random process, the RMS-slope, s, is related to the power spectrum 
by [19] 



s = {{C{x,)ff" = 5^p^^k^g{\k\), (26) 
and a good estimator for {D) has been shown to be [19] 



{D) 



!-^dkk-g{\k\) 



In the literature many different forms for the correlation function W(|a;i|) has been considered (see e.g. Ref. [6] and 
references therein). However, here we will only be dealing with two such forms. They are the Gaussian form given by 

W{\xr\) = exp(-^), (28a) 



9(1^1) = V^a exp ^- 



(28b) 
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FIG. 4: Examples of two rough profiles both with Gaussian height distributions and with an RMS-value & — 30nm. The 
power spectrum is of the (a) Gaussian type, with a — lOOnm, and the (b) West-O'Donnell type, with k- = 0.82(a;/c) and 
k+ = 1.97(w/c). Here the wavelength is A = 632. 8nm. With these parameters the RMS-slope and distance between consecutive 
peaks and valleys are respectively s = 0.424 and (D) = 128. 3nm for the Gaussian power spectrum, and s = 0.427 and 
(D) = 201. Inm in case of the West-O'Donnell power spectrum. Note that there are different scales on the first and second 
axes, with the result that the profiles appear much rougher than they are in reality. The two surface profiles where generated 
from the same underlying uncorrelated random numbers. 



where a is the transverse correlation length, and the so-called West-O'Donnell (or rectangular) form 

= «infc+xi-sinfc-:ri 

[k+ — k-)xi 

g(\k\) = —^!—[0{k+-k)e{k-k-) + e(k-+k)0(-k-k-)]. (29b) 
fe+ — k- 

where 9{-) is the Heaviside unit step function. In Eqs. (29) the quantities k±, with < A:_ < denote the lower 
and upper momentum cut-off for the spectrum, and they will be given a more precise definition in later sections. 
The latter power spectrum was recently used by West and O'Donnell [20] in an experimental study of the enhanced 
backscattering phenomenon from weakly rough surfaces. 

For the two above power spectra the mean slope, s, and the distance between consecutive peaks and valleys, (D), 
then become [19] 



s = 



{V^^, Gaussian 
73 



1 -^^k\ + k+k-+kl, West-O'Donnell ' 



and [19] 



^a, Gaussian 

{D) = { f^lJpF ^ „ • (31) 

, West-O'Donnell 

Two surface profiles with the same (Gaussian) height distribution, but with a Gaussian and a West-O'Donnell power 
spectrum possessing nearly the same value of the RMS-slope, s, are plotted in Figs. 4. 

It will later in explicate calculations prove useful to also have the Fourier representation of the surface profile 
function (and its inverse) at our disposal. They are defined as 

/°° tik 
^_C(fc)e^'=-S (32a) 
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C(fc) = / dxi C{xi)e-''''" . (32b) 

^ — oo 

The Fourier transform, C{k), of the surface profile function also constitutes a zero- mean Gaussian random process 
with statistical properties 

(C(fc)) = 0, (33a) 
(C(fc)C(fc')) = 27r6{k + k')6^g{\k\), (33b) 
where i5(-) denotes the Dirac delta function. 



2. Non-Gaussian Random Surfaces 

Naturally occurring surfaces have often more complicated height distributions then the Gaussian [6]. To fully 
characterize such surfaces are quite difficult and probably explains why they have gotten less attention in the literature 

then they probably deserve. The main problem is that in order to characterize them statistically, moments of in 
principle infinite order has to be known. These moments arc not, as for Gaussian surfaces, related to moments of 
lower order in a trivial way since the characteristic function is in general not known for non-Gaussian surfaces. We 
do not intend in this introduction to discuss non-Gaussian random surfaces in any detail, since we will not focus on 
them later. However, we would like to mention that as long as this kind of surfaces can be generated numerically, the 
scattering problem for non-Gaussian surfaces are not hard to handle by numerical simulations [21, 22]. On the other 
hand, small amplitude perturbation theory, say, can not be utilized in its standard form to non-Gaussian surfaces. 



G. Self-afflne surfaces 



It has been known for quite some time that self-afhne surfaces are abundant in nature. They can be found in 
various areas of natural science such as surface growth [23-25], fractured surfaces [26], geological structures [27, 28], 
metallurgy [29], and biological systems [30] to mention a few. 

A surface, Ci^i)-, is self-affine, according to its definition, between the scales ^_ and ^_|_, if it remains (either exactly 
or statistically) invariant in this region under transformations of the form 

Xi Xxi, (34a) 

C - A^C, (34b) 

for all positive real numbers A. Here H is the roughness exponent, also known as the Hurst exponent, and it charac- 
terizes this invariancc. It is usually found in the range from zero to one. When H = 1/2 the surface is an example 
of the famous random (Brownian) walk where the surface is uncorrected. However, if if ^ 1/2 the profile shows 
correlations; for H > 1/2 it is said to be persistent (correlated), and for H < 1/2 it is anti-persistent (anti-correlated). 
The reason for this naming is that if the self-affine "walker" when moving from the previous to the present space step 
went up, say, it is more likely that it will go up (down) in the next one ii H > 1/2 {H < 1/2). 
The scaling relation (34) is often put in the more compact, but equivalent form 

C(a;i) ~ A-^C(Aa;i), (35) 

where ~ is used to indicate statistically equality. This relation says that if we take the original profile C{xi)j enlarge 
(or contract) the lateral direction by rescaling xi into Aa;i, and simultaneously scaling to A~^C) the profile C{xi) 
and its rescaled version X~^({Xxi) should be indistinguishable. Of course, this holds true in an exact sense only for 
deterministic surfaces. In the statistical case, however, it is the statistical properties of the profile and its rescaled 
version that are indistinguishable. In Figs. 5 we presents some examples of self-afHne surfaces of Hurst exponent 
H = 0.3 (Fig. 5a), H = 0.5 (Fig. 5b), and H = 0.7 (Fig. 5c). As can be seen from these figures the landscapes become 
more "calm" the larger the Hurst exponent becomes. 

The scaling relation Eq. (35) does not fully specify the self-affine surface. In particular no information is contained 
in Eq. (35) about the amplitude of the surface. Such information is provided by the length scale, £, known as the 
topothesy. This length scale is define as the length, £, measured along the a;i-direction, for over which the root-mean- 
square of the height-difference between two points separated by £ is just £. To make this even more clear, let us 
introduce 



a(Axi) = ({C(xi+Axi)-C(xi)}')' , 



(36) 
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FIG. 5: Examples of self-affine profiles. The Hurst exponents were (a) H = 0.3, (b) H — 0.5, and (c) H = 0.7, and for all 
cases the topothesy where £ = 10~3A where A is an arbitrary length scale. The surfaces were generated by the Fourier filtering 
method from the same uncorrelated Gaussian distributed number. Notice how the RMS-height of the surface as measured over 
its total length increases as we increase the Hurst exponent. This is in agreement with Eq. (38). 



as the RMS- value of the height- difference measured over a window of size Axi . With this definition the topothesy is 
defined as the length scale for which 

a{e) = I. (37) 

From Eq. (36) it follows immediately that a{x\) ~ xf^, so that with Eq. (37) we get 

(T(Aa;) =^i-^Axf. (38) 

Notice that Eq. (37) allows for a geometrical interpretation of the topothesy as the length scale over which the profile 
has a mean slope of 45 degrees. The smaller I, the flatter the profile appears on a macroscopic scale. It should be 
stressed that in spite of the geometrical interpretation of ^, there is nothing a priori that restricts the topothesy to 
length scales where the self-affinity can be found. However, for the surfaces usually considered in scattering problems, 
we rather expect that f <C When ^_ < < ^_)_, the topothesy makes the transition between the scales, below 
for which a fractal dimension D = 2 — H can be measured using e.g. the box counting method [25, 31-33] and 
the scales, above £, for which this dimension is just unity. For length scales ^_ < Axi < £ the fractal dimension is 
therefore nontrivial (read different from one) and we have an example of a self-affine fractal [25, 31-33]. It should 
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be noticed that the fractal property of the self-afFine surface crucially depends on which length scale the surface is 
being observed. This essential point seems often to be overlooked in the literature where one too often treat self-affine 
surfaces as they were fractals [32, 33] at any length scale. 

Even if the self-afHne correlations of the profile is fully characterized by its Hurst exponent H, its topothesy 
parameter £ and the bounds of the self-affine regime ^_ and ^+ , nothing is said about its height-distribution. It is 
therefore not uncommon to talk about for example a Gaussian self-affine surfaces meaning that the surface correlation 
is of the self-affine type, while the distribution of heights is Gaussian. Thus by specifying the self-affine parameters, 
i.e. H, £, and in addition to the parameters needed in order to characterize the height-distribution, the surface is 
completely specified. Under the assumption that the surface has Gaussian height distribution it can be shown that 
the probability, p{C',Xi) for finding height ^ at position xi given that C(0) = 0, can be written as [25] 



1 



2Tr£ 



l-HrrH 



exp 



c 



(39) 



However, independent of the height- distribute being Gaussian or not, should satisfy the following scaling 

relation [25, 32] 



(40) 



which can be derived from the scaling relation Eq. (34). 

In fact, the scaling relation (34), or the equivalent form given in Eq. (35), is extremely powerful and can be used 

to derive most, if not all, of the properties of a self-affine surface. To show an explicit example of this, we would like, 
before closing this section, to derive the scaling relation of the power spectrum of the surface. This scaling relation is 
the most popular one to use for both generating self-affine surfaces as well as to measure the Hurst exponent. For a 
surface, C,{x\), of length, Li, the power spectrum is defined as 



5(1^1) 



= - T 



dxi e*'="MC(2/i+a;i)C(yi)) 



(41) 



where, as we recall, (C(|/i + a^i)C(2/i))yi is the (two-point) correlation function. By now taking advantage of the scaling 
relations (34) and (35), one finds 



1 



2 

Hence, one obtains from Eq. (42) that 



d(Axi) e''*-^(A^((yi+xi) A^C(yi)L ■ 



(42) 



(43) 



so that the power spectrum itself has to scale like 

g{\k\) ~ (44) 
For more details about self-affine surfaces and their properties the reader is referred to the literature [25, 32, 33] . 



H. Numerical Generation of Randomly Rough Surfaces 

Earlier in the sections we have discussed how to statistically characterize randomly rough surfaces. In analytical 
work, this is all what we need. However in a numerical Monte Garlo simulation approaches to be presented in a 
later section, individual surface, called realizations, have to be generated so that they possess the right statistical 
properties. The question therefore is: How can we do this? We do not intend to give a detailed discussion here, but 
will instead sketch how it can be done. 

As long as the power spectrum of the surface is known, an efficient way of generating the surface is by using the 
so-called Fourier filtering method [25, 32]. This method basically consists of to main steps. First, uncorrelaiedxm\Acra\ 
numbers of the type wanted for the height distribution of the surface are generated in Fourier space. Second, these 
numbers are filtered by the square root of the power spectrum 5'(|A;|), and the result transform by an inverse Fourier 
transform to real space. It was in this way that the surfaces shown in Figs. 4 and 5 were generated. For more details 
the reader is advised to consult Refs. [21, 22, 25] and [32]. 
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FIG. 6: The main scattering geometry used throughout this section for the wave scattering from a rough surface defined by 
^3 = C(3^i)- The region above the surface, X3 > CC^^i); is assumed to be vacuum {e+{uj) — 1), while the one below is metal or 
a dielectric characterized by a frequency-dependent dielectric function £-(0;) = £(aj). Notice for which direction the incident 
(^0), scattering {9s), and transmission {Ot) angles are defined positive. An angle of transmission is only well-defined if the lower 
medium is transparent, i.e. if Res{uj) > 0. 



III. QUANTITIES AND TECHNIQUES USED IN ROUGH SURFACE SCATTERING STUDIES 



The intention of the present section is to introduce some of the main quantities and techniques, both analytical and 
numerical, used in the field of wave scattering from randomly rough surfaces. The idea is to cover in some detail the 
most central techniques at the sacrifice of a wide coverage. The first part of this section is devoted to the discussion 
of some general properties of the scattering problem. This discussion is independent of the techniques used for its 
solution. In the second part, however, some of the central theoretical approaches towards the solution of the scattering 
problem are presented. Here the outlined theories will only sparsely be applied to a concrete problem. However, this 
is done in the next section where phenomena in rough surface scattering are discussed. 



A. The Scattering Geometry 



The scattering geometry that we will mainly concern ourselves with in this section is depicted in Fig. 6. It consists 

of vacuum (£+(0;) ~ Sq{lu) = 1) in the region X3 > C(-^i)j s,nd a metal or dielectric characterized by an isotropic, 
frequency-dependent, dielectric function S-(oj) = e{u>), in the region ^3 < Ci^i)- Here Ci^i) denotes the surface 
profile function and it is assumed to be a single-valued function of xi that is diflferentiable as many times as is 
necessary. Furthermore, it constitutes a zero mean, stationary, Gaussian random process which we from Sect. II F 
recall is defined by 

(C(xi)) = 0, (45a) 
(C(xi)C(a:;)) = S^W{\x^-x[\). (45b) 

Here W^(|a;i|) denotes the auto-correlation function and it will be specified later. 

The incident wave is assumed to be either p- or s-polarized, as indicated by the index v, and the plane of incidence 
will be the .x 1.9:3 -plane. Furthermore, the angle of incidence, reflection, and transmission, 60, 6s, and 6t respectively, 
are measured positive according to the convention indicated in Fig. 6. 



B. The Scattered Field 



From Sect. II we recall that in order to solve the scattering problem, we have to solve the Helmholtz equation and 
satisiy the boundary conditions, Eqs. (14), at the rough interface ({xi) as well as the boundary conditions at infinity. 
In the present section we will give the form of the far fields that automatically satisiy the Helmholtz equation and 
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the boundary conditions at infinity. We will discuss separately the case where the incident field is a plane wave and 
where it is a wave of finite width. 

However, before wo do so, wo rccali that for the scattering geometry depicted in Fig. 6, the Maxwell's equations 
are equivalent to the scalar Helmholtz equation for the field xslw) defined by Eq. (11), i.e. 

$.(.„X3|a;) = /^fi'-3N u = p, ^^g^ 

[ E2{Xi,X3\,UJ), V = S, 

It is the asymptotic, far-field behavior of ^v{xi,X:i\uj) that we are trying to determine. 

1. Plane Incident Wave 

Let us first consider the case where the incident field is (a either p- or s-polarized) plane wave of the form 

$t"°(a;i,X3|a;) = e*'=^^-*"«('=''^)^^ (47) 

where^ 



ao(g,w) = a+{q,uj) = { V } " . (48) 



Then the form of the field in region 2:3 > max^(xi) that satisfied both the Helmholtz equation as well as the boundary 
conditions at infinity (2:3 = 00) can be written as 



(49a) 



Similarly, a solution to the Helmholtz equation in the region X3 < m.m({xi) that satisfy the boundary condition at 
X3 = —00 is 

<p-{xuX3\u;) = ^ T,(g|fc)e^«^l-'«(«•'^)^^ (49b) 

where (e-(w) = s{w)) 



a;2 

a{q,u)) = a-{q,u)) = ye{u)——q'^, Rea,lma>0. (50) 

In these equations R^{q\k) and T^{q\k) denote the scattering and transmission amplitudes respectively. Notice that 
these asymptotic expressions does note say anything about how the fields look like in the surface region minC(a:;i) < 
X3 < maxC(a;i). This and its consequence, will be discussed in more detail in Sect. lllE when we derive the so-called 
reduced Rayleigh equation. 

2. Finite Width Incident Wave 

If the incident field is not a plane wave, but instead has a finite width, then the above expressions will have to be 
changed somewhat. In this case the incident field can be written as 

p(J.ykx,-^ao{k,u.)x3 ^ (51a) 



We will use the notation ao{q,u>) instead of a+{q,u>) in order follow the notation frequently used in the literature. 
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i.e. as a weighted sum of plane waves. Here F{k) is in principle an arbitrary function for which the integral exits. 
Due to the linearity of the Maxwell's equations, the scattered field becomes 



where 



f 

J —I 



dq 



R„{q,(j)e 



iqxi+iaQ{q,ui)x3 



\^R.{q\k)F{k). 



(51b) 



(51c) 



The total field in the region .T3 > maxC(.Ti) therefore is $+(xi, .T,j|Li,') = $™'^(a;i, a;3|a;) + $*'^(a;i,a;3|a;). 
In a similar way the field in the region 2:3 < minC(a;i) can be written as 



(52) 



where T,y{q,oj) is given by an expression similar to Eq. (51c). 

In order to fully define the asymptotic forms of the field, the envelope F(k) has to be given. Here we will only 
consider so-called Gaussian finite beams. Such beams are obtained if F{k) has the Gaussian form. If the half-width 
of the incident beam is denoted by w the Gaussian envelope F{k) can be written as [21] 



m = 



WU! 



2^/^^cao{k,u>) 



exp 



2 9 

W LO 



4c2 ( 



arcsm Oq 



(53) 



C. The mean differential reflection coefficient 



In the previous section, we obtained the asymptotic forms of the scattered and transmitted fields. These fields are 
known whenever the scattering and transmission amplitudes Rv{q\k) and Tv{q\k) are known. We will later in this 
section describe methods for how to determine these amplitudes. 

However, these two amplitudes are not accessible in experiments. Since, of course, our ultimate goal is to compare 
the theoretical predictions to those of experimental measurements, one has to relate these amplitudes to measurable 
quantities. Such quantities are provided by the so-called mean differential reflection and transmission coefficients. 
These are not the only experimentally accessible quantities possible. However, other such quantities must necessarily 
be related to the reflection or transmission amplitudes, since they fully specify the scattering and transmission problem. 

The mean differential reflection coefficient'' is defined as the fraction of the total incident power scattered, by the 
randomly rough surface, into an angular interval of width d6s about the scattering angle Og. Thus, in order to obtain 
an expression for this quantity one has to find an expression for the power incident onto the rough surface and the 
power scattered from it. We recall that the total power contained in an electromagnetic wave of electric and magnetic 
field vectors E and H respectively is given by the real part of the complex Poynting vector S = E x H*, where the 
asterisk denotes complex conjugate. More useful to us is in fact the time-averaged of this (complex) quantity. It is 
given by [10, 11, 34] 

(S), = iExH*, (54) 

where indicates time average. Hence the time-averaged power incident onto the rough surface, and scattered 
from it, are given by the real part of the 3-component of (S)j, evaluated for the fields involved. The corresponding 
time-averaged total energy fiux therefore becomes^ 

P = I dxxdx2 Re (S)^ = L2 / dxx Re . (55) 



If the surface is two-dimensional, which we however will not discuss in great detail here, one has to consider scattering into solid 
angle dSlg around the scattering direction {6s,<t>s) instead of into the angular interval d9s around the scattering angle 9s as is the case 
if the surface is one-dimensional. Furthermore, one also has to take into account that depolarizations may occur in scattering from 
two-dimensional surfaces. Hence the mean differential reflection and transmission coefficients in the 2D-case have polarization indices 
referring to the polarization of the incident and scattered light respectively. 

Recall that the coordinate system is chosen so that the a;iX2-plane coincide with the average (planar) surface, and thet the incident 
plane is the aJixs-plane. 
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In writing the above equation we have taken advantage of the fact that for a one-dimensional rough surface with its 
generator along the Xi-direction (as we consider here), the a;2-integration becomes trivial and only contributes with 
a factor Z/2, the length of the surface in the X2-direction. 



1. Plane incident wave 



We recall that if the incident wave is a plane wave of the form given in Eq. (47), then the scattered field is given 
by the second term of Eq. (49a). The incident and scattered energy fluxes thus becomes 

Pine = ^^^—o:o{k,Lo), (56a) 

Z CO 

and 

^^'^ ^ Tw" / „ If "o^*''^) l^^(9|fc)|% 

= I dOs Pedes), (56b) 

where 

Psc{Os) = ^iocos^ 9 s\Rv{q\k)\^ ■ (56c) 
47r 

Hence, the differential reflection coefficient, according to its definition, is given by the following expression 

dRy _ PscjOs) _ 1 uj cos^ 9s I p . I 

09 s ~ P^nc ~ il27rc COS ^^"^'^^"'^ ' 

In the above expression it is understood that the momenta k and q are related to the angles 9o and 9s according to 

k = — sm9o, (57a) 
c 

q = — sin 6*8. (57b) 

c 

Notice that dR,^/d9s includes the contribution from only one single realization of the rough surface. However, we 
are more interested in the mean of this quantity obtained by making an average over an ensemble of realizations 
of the rough surface profile. In consequence we obtain the following expression for the mean differential refl.ection 
coefficient (DRC) 

dRr,\ 1 U COS'^ 9s /,„ , ,,.|2\ i^^. 

-89:) - L^2;^W^(I^''(«I')I )• 

When light is scattered from a randomly rough surface both coherent (specular) and incoherent (diffuse) scattering 
processes will normally occur. The scattered power due to both these processes are contained in Eq. (58). In theoretical 

studies of wave scattering from rough surfaces it has proven useful to separate these two contributions even though 

such a separation is not possible under experimental conditions. The separation is done by noticing that ^|i?i,(q'|fc)| 
can trivially be written as 

R4q\k)\^) = {\RAq\k)\') - \{RMk))\' + \{RAQ\k))\^ • (59) 

Here the last term (on the right hand side) corresponds to the coherently scattered light, while the first two terms are 
related to th(^ light scattered incoherently. By using this result we find that the mean DRC can be subdivided into a 
coherent and an incoherent part, and they are respectively given by 

dR^\ 1 UJ cos2 6', 2 

\{R^{q\k))\ , (60a) 



|2 



d9s /coh -^1 2^'- cos^o 
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and 



dRij \ I to cos^ 9s 

/ incoh ^1 27rc cos 6*0 



RMk)\')-\{RMm' 



(60b) 



Expressions for the m,ean differential transmission coefficient, {dTi^/dOs) = {ptr{Ot) / Pine) , can be obtained in an 
analogous way by calculating ^^^(^'f), the equivalent oipsc{Qs) in transmission. The results of such a calculation is that 
the expressions for {dT^/dOg) can be obtained from those of {dRv/dOg) by substituting the transmission amplitude 
Ti,{q\k) for the reflection amplitude, Ri,{q\k) and multiplying the final expression with a factor of y^e(tJ). 



2. Finite width incident beam 



In the previous subsection we considered a plane incident wave. Such waves can never be achieved under experi- 
mental conditions, and it is therefore desirable in some cases to work with a incident beam of finite width. 

Such a beam has already been defined in Subsection. Ill B 2 and with these expressions one gets for the incident 
and scattered energy fluxes 

(61a) 

where erf (a;) is the error-function [39], and 

Psc{Os) = L2^cos^es\RM\k)\' ■ (61b) 
zvr 

Hence, the difi'erential reflection coefficient, according to its definition, is given by the expression 




(62) 

Also here q is understood to be related to the scattering angle by Eq. (576). The coherent and incoherent part of the 

mean differential reflection coefficient arc obtained in the same way as for the case of a plane incident wave. These 
expressions will not be explicitly included here since they follow from Eq. (62) by a simple substitution for of Eq. (59). 



wc 

L2 ;= 



D. General Properties of the Scattering Problem 

In order to solve the scattering problem, we have to calculate the scattering and transmission amplitudes, R,^{q\k) 
and T^{q\k). However, before we start discussing various methods for obtaining this goal, we will introduce some gen- 
eral features that the scattering problem should fulflll. These properties are, among others, reciprocity and unitarity. 



1. Reciprocity 

A general property of the scattering problem is the one of reciprocity. It involves the scattering matrix, Sj,((]\k), 
defined via the scattering amplitude according to 

SMk) = ^^pi?.(#). (63) 

The reciprocity theorem states that this scattering matrix, or just S-matrix for short, should satisiy the following 
relation 



Su{q\k) = S,{-k\-q). 



(64) 
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This relation for the surface scattering problem can under rather general assumptions be derived rigorously from 
Lorentz's reciprocity theorem [8]. However, we will not present such an interesting, but lengthy derivation here. 
We would, however, like to point out that such a derivation does not assume anything about the dielectric functions 
involved. Furthermore, there is no restriction on how strongly rough the surface is, neither how it is correlated. Hence, 
the reciprocity theorem is generally valid. It should also be pointed out that there seems to be no equivalent theorem 
to Eq. (64) that involves the transmission amplitude Tv{q\k). Reciprocity is therefore a property of the scattering 
amplitude. 



2. Unitarity 

In cases where the scattering medium is a perfect reflector, i.e. if Re£(w) < and Im£(ti)) = 0, the scattering 
matrix, <S'j/(g|fc), possess an additional property call unitarity. Since there is no absorption (Im£(a;) = 0) and no 
transmission in the system the energy incident on the rough, perfectly reflecting surface must be conserved. Without 
going into details, this has the consequence that the following relation has to be satisfied [8] 

/ 1 ^ = 27r(5(fc - k'), |fc|, \k'\ < ^. (65) 

It can be derived by calculating the total energy flux scattered from the surface that, due to energy conservation, 
should equal the incident energy flux. Eq. (65) expresses the unitarity of the scattering matrix, and it is a consequence 
of the conservation of energy in the scattering process. 

Even if Eq. (65) is derived under the assumption that energy conservation is satisfied, lot us for a moment show 
how this indeed follows from the unitarity condition. Let us assume that the rough surface has length Li. Then 
we know from the sampling theorem [40] that the smallest momentum variable that we can resolve is 2'k/Li. By 
multiplying each side of Eq. (65) by rffc'/(27r) and integrating the resulting expression over an interval of length 2'k/Li 
that contains k' = k one finds 

hFM = ^- ^^^^ 

By now using the definition of the scattering matrix, Eq. (63), together with Eq. (576), one arrives at 

' ^^'^ \ d9s = 1, (67) 



n/2 \ d9s 

where we have taken advantage of Eq. (58) . From the definition of the mean differential reflection coefficient given in 
Sect. IIIC 1, we understand that Eq. (67) is just the conservation of energy for the scattering system considered. 



3. Energy Conservation 

If, however, the lower medium is not a perfect conductor, but still is a non-absorbing medium (Im£(a;) = 0) the 
unitarity condition Eq. (65) will no longer hold true. However, we should still have conservation of energy. This 
means that all energy incident onto the rough surface should be either scattered from it or transmitted trough it. 
This fact is expressed by the following equation 

u!,'={eo,oj)+ul'{eo,oj) = i, (68) 

where is the angle of incidence of the light, and 



it/2 



OR, 



K'{Oo,co) = I ^^{^-^) de,, (69a) 
U'Jieo,^) = 1^) dOt. (69b) 



J —IT 



/2 \ det 

Physically Uf,'^{6o,uj) expresses the fraction of the incident energy scattered from the surface, while in a similar way 
UI'{Oq,uj) expresses the energy fraction transmitted through the system. Notice that the energy conservation condition 



23 



should hold true for all incident angles and polarizations as well as being independent of the width of the incident 
beam. The only restriction being than there should be no medium that absorbs energy. However, if absorption 
is present, Eq. (68) might be modified by adding an absorption term to the right hand side. Unfortunately, this 
absorption term is hard to calculate in a rigorous way. 

For practical purposes, the conditions Eqs. (67) and (68) are most frequently used as a test of the quality of numerical 
simulations (see Sect. Ill I). In such an approach these conditions are necessary, but not sufficient conditions for the 
correctness of the simulations. 



E. Derivation of the Reduced Rayleigh Equation 

The reduced Rayleigh equation (RRE), under which name we know it today, was first derived by Toigo, Marvin 
and Celli [44] in the last half of the 1970's. This equation is the single integral equation satisfied by the reflection or 
transmission amplitudes. This equation, even if its precise region of validity is hard to quantify in detail, has served 
as the starting point for many, if not all, of the perturbative techniques developed in the field of wave scattering from 
rough surface. We would, however, already at this early stage like to stress that the reduced Rayleigh equation is not 
restricted to the same limitations as perturbation theory, and that its validity goes beyond that of such theories. It 
can in fact also be used in numerical simulations to obtain non-perturbative results. 

We will below give the detailed derivation of the RRE for reflection. The scattering geometry that we consider 
is the one presented in Fig. 6. This geometry is illuminated from above by a plane incident wave of either p- or 
s-polarization, and the incident plane is the xiccs-plane. 

1. The Rayleigh Hypothesis 

It should be clear that for the region above the maximum point of the surface the total field takes the form given by 
Eq. (49a) and similarly the total field below the minimum point of the surface can be expressed according to Eq. (496). 
However, in order to solve the scattering problem, one has to take into account the boundary conditions to be satisfied 
at the randomly rough surface C{xi) separating the two media above and below it. The problem is that we don't know 
the form of the total field close to the surface, or to be more precise, in the region min^(xi) < ^(2^1) < max^(xi). 
It should be obvious from a ray optical point of view, that at least for rather rough surfaces, expansions of the form 
(49) are not adequate to describe the total field in this region due to the lack of not allowing downward propagating 
scattered modes. However, as the surface becomes smother and smother, the asymptotic expansions of the field given 
earlier should represent a better and better approximation for the total field. This lead Lord Rayleigh [1, 2], when 
studying scattering from sinusoidal surfaces, to assume that the asymptotic expansions for the total field was not only 
valid in the region far away from the rough surface, but could also be used all the way down to the rough surface. 
Under this assumption, known today as the Rayleigh hypothesis, he could satisfy the boundary conditions on the 
rough surface and thereby derive equation which lead to the solution of his scattering problem. 

The Validity of the Rayleigh Hypothesis 

We will in the next subsection demonstrate this procedure when applied to the wave scattering from a randomly 
rough surface. However, before we do so, we will dwell a little upon the validity of the Rayleigh hypothesis. Theories 
based on this approximation do not properly include, as mcntionc;d above, downward propagating scattc;rc;d or upward 
propagating transmitted waves. From a naive geometrical optics argument, we realize that a scattering process 
producing incoming scattered or transmitted waves has to be a multiple scattering process. So, scattering geometries 
where the Rayleigh hypothesis is not valid, has thus to be dominated by multiple scattering, and has therefore to be 
rather rough. It should, however, be stressed that this do not imply that for any scattering geometry dominated by 
multiple scattering, the Rayleigh hypothesis is doomed to break down. One might very well have processes dominated 
by multiple scattering without receiving essential contributions from downward propagating scattered waves. Good 
examples of this arc provided by the ability of perturbative and numerical studies based on the Rayleigh hypothesis 
to show multiple scattering phenomena like the enhanced backscattering and satellite peaks [21, 42, 121]. 

The Rayleigh hypothesis is hence a good approximation if the surface is not too rough. H owever, at what level of 
roughness must we say that this approximation no longer is valid? There has been many papers devoted to the study 
of the validity of the Rayleigh approximation. There seems today to be consensus on the criterion [43] 



(70) 
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where S and a are the RMS-height and correlation length of the surface respectively. The reader is encouraged to 
consult the literature for more details [43]. 



2. The Rayleigh Equations 



In this subsection we will derive a set of two inhomogeneous coupled integral equation for the scattering and 
transmission amplitudes, Ri,{q\k) and Ti,{q\k). These equations are referred to as the Rayleigh equations, and we will 
now demonstrate how they are obtained. 

From Sect. II C we recall that the boundary conditions to be satisfied by the field on the surface are the continuity 
of the field and its normal derivative, i.e. 



dn^u (a;i,a;3|a;) 



X3 = ({xi) 



(71a) 
(71b) 



where the normal derivative, dn, and the symbol k^(u;) have been defined earlier in Eqs. (146) and (lAd). 

If now the Rayleigh hypothesis holds true, the asymptotic field expansions of Sect. IIIB, can be used in order to 
fulfill those boundary conditions. By substituting the asymptotic expansions, Eqs. (49), into the boundary condition 
for the field, Eq. (71a), one is lead to the following integral equation 



Akxi — iao{k 



iqxi — ia{q, u})C{xi) 



If we now rewrite this equation, which will prove useful later on, by using the properties of the Dirac 5-function we 
get 



^ „iqxi 



2n 



2Tr6{q - fc)e-*"°(9'")^(^i) + R4q\k)e'°'° 



dq 
2^ 



e''"'' T^{q\k) e-*"(«''^)';(=^i). 



(72a) 



By doing the same for the boundary condition for the normal derivative, i.e. Eq.(716), one arrives at 
' dq 



fdq 
J 2tt 



^iqxi 



-2TTS{q - k) {C'{xi)q + ao(g, ^)} e-'"''(9''^)«^i) + R,{q\k) {-C,'{xi)q + ao((Z, w)} e^^^^'i'^^^^'''^ (72b) 



-ia(g,aj){(a;i) 



Together Eqs. (72) constitute a set of coupled inhomogeneous integral equations announced earlier — the Rayleigh 
equations. 



3. The Reduced Rayleigh Equations 

We will now continue to derive the so-called reduced Rayleigh equation (RRE) for reflection [44] and transmis- 
sion [45]. The RRE is a single integral equation satisfied by the reflection or transmission amplitude. They are 
derived from the Rayleigh equations by eliminating respectively the transmission and reflection amplitudes. 



The Reduced Rayleigh Equation for Reflection 



In order to obtain the reduced Rayleigh equation for reflection, we have to eliminate the transmission amplitude, 
Ti,{q\k) from the (coupled) Rayleigh equations given in the previous subsection. By multiply Eq.(72a) by 

g-ipxi-ia(p,a.)C(xi) [-Q'i^xx)p + a{p, w)] , (73) 

and Eq.(726) by 



(74) 
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adding the two resulting equations, and integrating the final result over .xi, one finds that the terms containing the 
transmission amplitude vanishes identically. In detail what one gets for terms proportional to T^{q\k) are 

[dxi^T^{q\k) [-Cixi){p + q) + a{p,uj)-a{q,ij)] e-»(p-9)^ie*(-«(f'^)-«(«'^»^(^i). (75) 
J 27r 

This expression is simplified by introducing an integral defined according to^ 

7(7|g) = y"da;i e-*T^(^i)e-*«^\ (76) 



From this definition it follows that 

g-f(7|g) 
7 

With Eqs. (76) and (77), Eq. (75) can be written in the form 

{p- q){p + q) 



= y^da;i C'(a;i)e-'Tf(^i)e-*«^\ (77) 



a{p,u>)+a{q,(j) 



+ a{p,Lo) -a{q,u>) 



I{a{p,uj) + a{q,ui)\p - q). 



By some simple algebra it can readily be shown that the expression in the square brackets is identically zero, and 
thus the transmission amplitude T^{q\k) has been eliminated from the Rayleigh equations. 

The reduced Rayleigh equation for reflection now follows from the remaining non-vanishing parts of the equation. 
It reads 

J g M+{p\q) RMk) = M-{p\q) (78a) 



where 

' {p + Kt,{uj)q){jp - q) 



Mt{p\q) = ± 



a{p,u)Tao{q,bj) 



+ a(p, Lo) ± Ky{uj)aQ{q, uj) 



P-polarization 



J(a(j3,w) T ao(g,'^)b- g) • (78b) 



If we restrict ourselves to p-polarization, the reduced Rayleigh equation as presented in Eqs. (78) takes on a simpler 
form. By a straight forward calculation one finds^° 

J^N+ip\q)Rp{q\k) = N-{p\q), (79a) 



where 



Np{p\q) = ± — -, r— -, r-I{a{p,u})Tao{q,u))\p-q). 

a[p,uj)Tao[q,uj) 



S-polarization 



(79b) 



Similarly for s-polarization one finds^^ 

j^Nt{p\q)Rs{q\k) = N-{p\q), (80a) 



where 



N^ip\q) = ±- ^ -. -Iia{p,uj)Tao{q,u;)\p-q). (80b) 



^ Be aware that various sign conventions seems to appear in the literature for this quantity. 

Here the matrix elements M^{p\q) and N^{p\q) are related by M^(p\q) = (e(a)) — l)N^{p\q). 
11 Here is M±(p|g) = {e{uj) - l){uj/c)''N^(j,\q). 
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4- The Reduced Rayleigh Equation for Transmission 

In the previous subsection it was shown that by eUminating the terms containing the transmission amphtude from 
the two coupled Rayleigh equations, the reduced Rayleigh equation for reflection was obtained. In a similar way. the 
scattering amplitude might be eliminated from the same equations, resulting in the reduced Rayleigh equation for 
transmission. 

We will not go into details about how this is done, since the derivation mimics the one given in the previous 
subsection. Here we will only give the results that are [45] 

f°° dq I{a{q,uj) ~ aa{p,u;)\p- q) M^/ i/x o , .2e{u>)aQ{k,oj) 

/ 7i 7 ^ 7 ^ m + <^o{p,oJ)a{q,uj)]Tp{q\k) = 2-k5{p - k) \\ \ ■ (81a) 

for p-polarization, and [45] 



J — C 



dqI{a{q,i.)-a,{pM\p-<l) ^^^^^^^ ^ 27r%-fc)- ^""^'^'^^ 



,27r a{q,uj) - ao{p,oj) '^{e{uj)-l)' 

(81b) 



for s-polarization. 

This concludes the section on the reduced Rayleigh equation. 



F. Small Amplitude Perturbation Theory 

Among the oldest theories addressing rough surface scattering we find the small amplitude perturbation theory [4] . 
The starting point for this perturbation theory, like most of the perturbation theories developed for handling wave 
scattering from rough surface, is the reduced Rayleigh equation (78). If the rough surface is weakly rough, most of the 
light incident upon it is scattered into the specular direction. However, due to the surface roughness, a small fraction 
of the incident power is scattered away from the specular direction. Theoretically, this non-specular scattering is 
taken into account by assuming an expansion for the scattering amplitude in powers of the surface profile function of 
the form^^ 

i?.(#) = f:^^. (82) 

n=0 

Here R^\q\k) is assumed to be of order 0(C") in the surface profile function Q{xi). In order to solve the scattering 
problem in this way, we therefore have to determine the expansion coefficients {R\'^\q\k)\. However, to determine 
all these coefficients is obviously not practically possible if C,{xi) is a rough surface. Therefore, the expansion (82) 
is terminated at some upper value A'', resulting in an A'''th order perturbation theory. In practical application one 
usually has N < 3-5. If the surface is weakly rough the sum of these N terms will provide a good approximation to 
the total scattering amplitude R,y{q\k). However, as the surface roughness becomes stronger and stronger, a higher 
number of terms has to be included in the expansion, and the method becomes cumbersome and not very practical 
since R^iT^qlk) for big values of n easily becomes complicated. Thus, the small amplitude perturbation theory is 
only of interest for weakly rough surfaces. Hence, it should therefore not represent any restriction to assume that the 
Rayleigh hypothesis is valid and therefore that R^{q\k) should satisfy the reduced Rayleigh equation, Eq. (78). 

Under this assumption, the various terms in the expansion for the scattering amplitude can in principle be obtained 
by substituting its expansion, Eq. (82), into the reduced Rayleigh equation, Eq. (78), and satisfy the resulting equation 
order- by-order in the surface profile function Ci^i)- However, before this can be done, N^{p\q), or equivalently 
M^{p\q), that enter via the reduced Rayleigh equation, also have to be expanded in powers of the surface profile 
function. Since the matrix-elements, N^{p\q), only depend on ({xi) through the integrals I{j\q), defined in Eq. (76), 
one makes the following expansion 

I{j\q) = /dxi e-*^'^("i)e-'«"i = C'"^?)^ (83a) 

71=0 ^' 



The prefactors in this expansion is included for later convenience. 



27 



where 

/oo 
dxiC{xi)e-"'''\ (83b) 
-oo 

is the (inverse) Fourier transform of the nth power of the surface profile function. From the above equations it should 
be apparent why we made the choice we did for the prefactors in the expansion for R^{q\k). 

When the expansion (83) is substituted into the expressions for N^{p\q), Eqs. (796) and (806), one obtains 



n=0 

where for p-polarization 



(84a) 



and 

Ar±(")(p|g) = {-ir[±pq + a{p,cu)ao{q,iv)][a{p, CO) Taoiq,u;)r-\ (84c) 
when n > 1, while for s-polarization we have 

and 

A/;±(")(p|5) = ±{-ir[a{p,u;)Tao{q,u;)r-\ (84e) 

when n > 1. 

With these relations available, a recurrence relations for {Rl^\p\q)}, is readily obtained by substituting Eqs. (82) 
and (84) into the reduced Rayleigh equation, Eqs. (79a) and (80a), and equating terms of the same order in ({q). 
The recurrence relation reads 

El") /"?-^.+ ^"""H'?b)C'"-"^(<?-p)4'"Hp|A;) = AC^"Hal^X^"H'/-fc)- (85) 

Now the expansion coefficients for the scattering amplitude, R^\p\q), should be rather straight forward to obtain 
(at least in principle). The lowest order term, n = 0, is given by 

Rl°\q\k) = 2Tr 6{p-k)R^°\k,oj) (86a) 

where R^\k,uj) given by (e = 1) 

' £{oj)aa{k,uj) - eo{uj)a{k,oj) ^ ^ 
s{u>)ao{k,Lu) + soito)a{k,u))' ' /ggj^x 

ao{k,uj) - a{k,uj) ^ > 

^ ao{k, w) + a{k, uj) ' 



In these expressions, the momentum variables k and q are understood to be related to the angles of incidence and 
scattering, 6o and 6s respectively, through 

k = — sin^oi (87a) 
c 

q = -smOs. (87b) 
c 

The above results are, as expected, the result that we would get for the scattering from a planar surface, and we note 
that R''u\k,uj) is nothing else then the Fresncl reflection coefficients [10]. Notice that the (5-function in Eq. (86a), 
coming from (^(°) , guarantees that the scattering is only into the specular direction. 
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The results for the higher order terms (n > 1) in the expansion of R^(q\k), describe the hght scattered by the 
roughness into directions other than the specular. These terms can be calculated recursively from Eq. (85), but, 
unfortunately, such expressions easily become rater cumbersome for higher the order terms. However, the first few 
terms are manageable, and they are given by the following expressions [46] 



2! 



xi'\q\p\k)C{q-p)ap-k), 



(88a) 
(88b) 



R^^\q\k) = ^ /_°° ^ ^X^^\q\Pi\P2\k) C{q - Pi)C{Pi - P2K{P2 - k), (88c) 

where, the functions Xi^\q\k), X!?\q\p\k) , ■ ■ ■ are somewhat lengthy functions of their arguments, and arc therefore 
given separately in Appendix B. 

As discussed earlier, the first few terms of the expansion (82) with Eqs. (88) substituted, should hence for a weakly 
rough surfaces represent a good (3rd order) approximation to R„{q\k). Experimental accessible quantities could 
therefore be calculated based on this approximation. For instance, the mean differential reflection coefficient from the; 
incoherent component of the scattered field is then, according to Eq. (606), through fourth order in the surface profile 
function, given by 



^ I incoh 



1 uj cos Bs 

L\ 2-KC COS^o 



i?«(5|fc) 



+ 



+ ii?e(i?i,3)(g|fc)i?W*(g|fc) 



i?(2)(g|fc)^ 

+ 0((5^). 



i?(2)(g|fc) 



(89a) 



In arriving at Eq. (89a) it has been assumed that the surface profile function, C,{xi), constitutes a zero- mean, stationary, 
Gaussian random process. Due to the Gaussian character of the surface only terms that contain an even number 
of surface profile function survives the averaging process. The different averaged contained in Eq. (89a) can all be 
related to the set of functions {)(^\q\k)} according to 



i?(2)(9|fc) 



i?(i)(<z|fc) 

R'^^\q\k) 



L,6^g{\q^k\) x«(«|fc) 
dp 



(89b) 



= Li5^ 



■.2ti 



9{\q - P\)9{\P - k\) 



XS'^ {q\p\k) + xS'^ («b|fc)xL'^ *{q\q + k-p\k)\, 



(89c) 



and 



{l6^\q\k)R^^^*{q\k)) = LiS^g{\q-k\)x';}^*{q\k) 

^ {xl;'\q\p\q\k)9{\p - q\) + X^^Hq\k\p\k)g{\p - k\) + x^^Hq\p\p + k - q\k)g{\p - q\)} . (89d) 

Notice that the term (896) represents single scattering, while the terms in Eq. (89c) give the contribution due to 
double scattering. Eq. {89d) represents a "mixed" contribution. 



Accuracy of the Small Amplitude Perturbation Theory 



So what accuracy can we expect to achieve by using the small amplitude perturbation theory, and for what range 
of surface parameters is it valid? There have been many studies in the past, both theoretical, numerical, and 
experimental, regarding this issue. Probably the most relevant for this introduction is the one by E. I. Thorsos and D. 
R. Jackson [57] who for an acoustic scattering problem studied the validity of small amplitude perturbation theory for 
Gaussian surfaces by comparing its prediction to the one obtained by rigorous numerical simulations (see Sect. Ill I). 
They found that the small amplitude perturbation theory is good if A;(5 <C 1 and ka < 1, where 6 and a respectively 
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arc the RMS-height and the correlation length of the surface. For an electromagnetic scattering problem this translate 
into 

^/\e{oj)\'^S < 1, (90a) 

and 

-a < 1, (90b) 
c 

The first condition (90a) comes from the fact that quantities containing the surface profile function should be expand- 
able in a Taylor series about their mean surface. The second criterion stating that the correlation length of the surface 
can not be too big originates from the fact that if the correlation length of the surface is too largo the Gaussian height 
distribution becomes close to a (5-function with the result that the second order term in the perturbative expansion 
can be of the same order as the first order term even if Eq. (90a) is satisfied. 

If the surface is non-Gaussian it is found that the above criteria can be relaxed somewhat [58] 



G. Unitary and Reciprocal Expansions 

In the previous section we presented small-amplitude perturbation theory, which is based on the expansion of 
the scattering amplitude Ri,{q\k) in powers of the surface profile function C{xi). In Sect. IIIDl we claimed that 
a valid theory for rough surface scattering should satisfy reciprocity, i.e. the theory should satisfy the relation 
Si,{q\k) = Si,{—k\ — q), where Su{q\k) is the scattering matrix as defined in Eq. (63). By inspecting the formulae 
obtained in the previous section for Rv{(i\k), it is at least not obvious that reciprocity is satisfied. Does this mean 
that small amplitude perturbation theory does not respect the principle of reciprocity, and therefore is an incorrect 
theory? The answer to this question is no. as you might have guest, and the small amplitude perturbation theory docs 
in fact respect reciprocity. However, to see this is not straight forward since an extensive rewriting of the expressions 
are need for. Theories where reciprocity is not apparent at first glance is normally referred to by saying that the 
theory is not manifestly reciprocal. 

Due to the lack of manifest reciprocity in the small amplitude perturbation theory as well as the desire to map the 
classical scattering problem onto the formalism of a quantum mechanical scattering problem [47], Brown et al. [48, 49], 
in the fist half of the 1980's, constructed a theory which was manifestly reciprocal. This theory goes today under 
the name of many-body perturbation theory, but is also known as self-energy perturbation theory. It is this kind of 
perturbation theory that we will concern ourselves with in this section. 



1. The Transition Matrix 



The starting point for the many-body perturbation theory is to make the postulation that the scattering amplitude 
Rv{q\k) should satisfy the relation [49] 

RM\k) = 27Td{q-k)Rl°\k,uj)-2iG^°Hq,uj)%{q\k)G^°\k,ij)ao{k,ij). (91) 

Here R^\k,(jj) is the Fresnel reflection coeflacient and defined by Eq. (86&). The second term of this equation, 
containing the transition matrix %{q\k), also known as the T-matrix, represents the field scattered away from the 

specular direction. Furthermore, Gl^\k,Lu) is the surface plasmon polariton Green's function for the planar vacuum- 
metal interface. This Green's function can be defined from the relation 

2iao{k,u)G^°\k,oj) + Rl°\k,oj) = -1, (92a) 

That leads to the following expressions 

is{uj) _ 

e(w)ao(fc,w) +a(fc,w)' ~ ^' (92^) 
aoikyto) + a{k,u)) ' v — s. 
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2. The Scattering Potential 
The transition matrix is postulated to satisfy the following equation [49] 

%(p\k) = V.{p\k) + p ^V.{p\q)G^^\q,uj)%{q\kl (93a) 
= VAp\k)+ f ^%{p\q)G^^\qMVu{q\k), (93b) 

where V{q\k) is known as the scattering potential. It is supposed to be a non-resonant function of its arguments, 

i.e. not containing the Green's function G^^\ In arriving at Eq. (936) we have used explicitly that both Vi,{p\k) and 
Ti,{p\k) are reciprocal, i.e. that Vi,{p\k) = Vu{—k\ — p) with a similar expression for the transition amplitude-^^. 

We now seek an (integral) equation satisfied by the T-matrix. This is done by substituting Eq. (91) into the reduced 
Rayleigh equation (78), and thereby obtain 

d<Z . I .^0. ^rr(\l.^ R^\k,oj)N+{q\k)-N-{p\k) 

2iao{k,uj)G^^\k,uj) 

Even though, the many-body perturbation theory could have proceeded from this equation, it has proven useful from 
a purely algebraic point of view, to instead of the T-matrix work in terms of the scattering potential Vv{(i\k). Thus 
we aim to obtain an integral equation for this quantity which our perturbation theory will be based directly upon. 
By substituting the right hand side of Eq. (93a) into Eq. (94), making a change of variable and using Eq. (93a) once 
more, one obtains the desired integral equation for the scattering potential. It reads 



/ 



^ A,{p\q)VMk) = B,{p\k), (95a) 



where the matrix elements are given by 



2iao{q,oj)G^°\q,u;) + R^J'\q,uj) N+{p\q) - {p\q) N+(p\q) + N-(v\q) 



2iao{k,Lo) 2iao{q,u)) 
and 

= R^\p,^)N-ip\k)-N,ip\k)_ 
2iao{k,uj)G)°\k,to) 

To obtain A,y{p\q) we have explicitly taken advantage of Eq. (92a). With the expression presented for the matrix 
elements for the reduced Rayleigh equation, Ni^{q\k), it is now straightforward to obtain closed form expressions for 
A^[p\q) and B^{p\k), but we will not present such expressions here. 

The integral equation (95) will be the starting point for our manifestly reciprocal many-body perturbation theory. 
The essence of the theory is to expanding the scattering potential in powers of the surface roughness C{xi ) according 
to 

n=0 

with similar expansions for Ai,{p\q) and B^{p\q). In these expressions the superscripts denotes, as earlier, the order 
of the corresponding terms in the surface profile function. 

One of the advantages of this theory is that even a lower order approximation to the scattering potential corresponds 
to a resummation of an infinite number of terms in an expansion in powers of the surface profile function^^ 



That this is indeed the case might be confirmed from the expressions to be derived later for tliesc quantities. 
This is also the case for the so-called self-energy perturbation theory to be presented in the next section. 
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We will not go into details here, but it can be shown that the results for the few first terms in the expansion of the 
scattering potential are [55] 



Vy\q\k) 



e{Lu)qk — a{q, u;)a{k,u))] Q^\q — k), v = p 



z^(£(a;)-l)CW(g-fc), 



V = s, 



(97a) 



for the 1st oder term, and 
V^''\q\k) 

and 



,. gH-i 

.(£H-1)2 



+2i 



[a{q, uj) + a{k, w)] [qk - a{q, u>)a{k, w)] C^^^ {q - k) 

a{q,oj) 



^C(i) (g - p)a{p, u;)C^'^ {p - k)a{k, oj), 



Vj'Hqlk) 



{s{uj) - 1) [a{q,uj) + a{k,w)] C^^\q - k), 



(97b) 



(97c) 



for the second order terms. Higher order terms can be found in e.g. Ref. [51]. It should be mentioned that by 
definition the lowest non-vanishing order of the scattering potential is 1st order in the surface profile function. In 

other words vj^\p\q) = always. 

As the reader easily may check, the above expressions are manifest reciprocal, i.e. 

V(/\q\k) = Vj^H~k\-q), 

and this property should hold true to all orders in the surface profile function [49]. However, it should be emphasized 
that expressions of this form are not obtained directly from the solution of Eq. (95), but that some rewritings instead 
are needed for [55]. 

If the transmission matrix is expanded in the same way as in Eqs. (96), i.e. 



%{q\k) = 5^i_LTj")(g|fc), 



(98) 



n=0 



a recurrence relation for {TS^\p\q)} in terms of Vu"^'{p\q) can be derived and thus T^{q\k) to some order can be 
calculated. 

Hence, through the calculation of the scattering potential the T-matrix is now known. The contribution to the 
mean differential reflection coefficient from the incoherent component of the scattered light is 



dR,. 



d0s 



^ I incoh 



1 2 /u;\3 

il TT V C / 



cos^ Os cos ^0 



q,LO] 



T.(g|fc)|^\ - |(T.(9|fc))|'l G^^\kM 



(99) 



This expression is obtained by substituting Eq. (91) into the defining expression for the incoherent component of the 
mean DRC (Eq. (606)). The expression in the square rackets for the 



H. Many-Body Perturbation Theory 

The Green's function Go((7,aj), we recall, is the surface plasmon polariton Green's function at a planar interface. 
In addition to this Green's fmiction it is also useful to define a rough surface Green's function, G,y{k,u)), or more 
formally the Green's function of a i^-polarized electromagnetic field at the randomly rough interface. Some authors 
refer to this function as the renormalized Green's function. It is defined as the solution of the following equation [49] 

G,{q\k) = 27r6{q-k)Gl°\k)+G^^\q) ^ ^VMp)G.{p\k). (100a) 

This equation is often in the literature referred to as the Lippmann-Schwinger equation for the renormalized Green's 
function G^,{q\k). Notice that from Eqs. (91) and (100a) it follows that R{q\k) = -2T:5{q - k) - 2iG{q\k)ao{k) . An 
alternative way of expressing the above equation is obtained by iterating on G^(g|fc). The result can be written as 



GMk) = 2n6{q-k)G<;°\k) + G<;^\q)%{q\p)Gl^\k), 



(100b) 
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where we have used a Born scries [47] expansion for the of the T-matrix in Eq. (93). This equation is often for 
simphcity expressied in operator form hke G = Gq + GqTGq. 

In terms of the renormaUzed Green's function the mean DRC takes on the form 



(^) = y-^S ^o^'^^^os^o \(\G.mf) - liGAiim'] , (101) 

/incohr Lnv L\ / J 



where the method of smoothing [50] has been appUed to Eq. (91) as well as the reduced Rayleigh equation. In these 
expressions the mean of the renormalized Greens function satisfies the Dyson equation 

{G,iq\k)) = 277,5(5 - fc)GW(fc) + C?W(g) J°° ^ (M.(g|p)) (G^k)) , (102a) 
where the unaveraged proper self-energy Mi,{q\k) is a solution of the equation 

M4q\k) = V.{q\k) + ^ M.(g|p)G(")(p) [V^k) - {V.{p\k)}] . (102b) 

Since the surface profile function is stationary, both the renormalized Green's function and the proper self-energy are 
diagonal in the momentum variables q and fc, i.e. 

{G,{q\k)) = 27TS{q-k)G4k), (103) 

with a similar expression for the proper self-energy. Under this assumption the renormalized Green's function can 
formally be written as 

= • (104) 

{G^^\k))-^ - M^k) 

Hence the surface polariton poles in Gi,{k) are shifted as compared to those of gI^'' due to the presence M^{k). The 
self-energy can be calculated perturbatively as an expansion in powers of the surface profile function. The resulting 
perturbation theory is known as self-energy perturbation theory [52] . 

The two-particle average Green's satisfies a Bethe-Salpeter equation [47] of the form^^ 

{\G.iq\k)f) = L,2n6iq-k)\G.{q)\' + \G.{q)\' ^U,iq\p)(\GMk)f) , (105) 
where U,y{q\p) is the so-called irreducible vertex function. Formally, one may write the solution to Eq. (105) as 
(]G,(g|fc)|') = L^27rS{q- k)\G^{q)f + L,\G^{q)f X4q\k)\G^{k)f 

(106) 

where Xjj{q\k) is the reducible vertex function. With this equation the mean DRC takes on the form 

^) = -^cosHsCOseo\G4q)fX4q\k)\G4k)\\ (107) 

'^''s / incohr ^ 

The reducible vertex function can be shown to be related to Ui,{Q\k) though the equation 

XM^) = Milk) + f ^ U4q\p) \G4p)f X,{p\k). (108) 

Unfortunately we do not know, in general, how to solve the Bethe-salpeter equation (105). Hence some approximative 
methods have to be employed. The most frequently used methods are the Preilikher factorization [53] , or a diagram- 
matic method [47, 54, 73]. In this latter approach X,,{q\k) is approximated by a subset of (an infinite number of) 



In arriving at this equation also here the method of smoothing [50] has been apphed. 



33 



diagrams. Those usually arc the ladder diagrams and the maximally-crossed diagrams, where the former describes 
wave diffusion in the random media, while the latter is related to wave localization. 

With this approach it can be shown that the reducible vertex function for p-polarized light incident on a rough 
metal surface can be written as [54] 



A 



( q-k I k-q \ 
\ 2 \ 2 ) 



4A2 (g + A;)2+4A2 



(109) 



which when substituted into Eq. (107) defines the mean DRC in this approximation. Here, A{q\k) is a smooth function 
of its arguments and A = + Agp is the (total) decay rate of surface plasmon polaritons due to Ohmic losses (A^) in 
the metal and conversion into other surface plasmon polaritons (A^p ~ Im{M{ksp))). Their mathematical expressions, 
as well as the other quantities appearing in Eq. (109), can be found in Ref. [54]. The first term in Eq. (109) is due to 
single-scattering, the second arises from the ladder diagrams, while the last one is the contribution from the maximally- 
crossed diagrams. This last term is the one that is responsible for the enhanced backscattering phenomenon that we 
will discuss in the next section. 

It should be mentioned that in arriving at Eq. (109) the pole- approximation for the renormalized Green's function 
has been utilized. This approximation amount to writing 



C(u;) 



■ ksp + iA 



(110) 



where ksp is the wave vector of the surface plasmon polariton (See. Eq. (19)), and C{u)) is a constant. The Green's 
function for s-polarization does not have poles, and the pole-approximation is thus not relevant in this case. 



I. Numerical Simulation Approach 

In the previous sections various perturbation theories were discussed. Such theories catch the main physics of 
the scattering problem if the surface is not too rough. However, for interfaces that are strongly rough, none of the 
perturbative approaches can be used because too many terms in the expansion have to be included in order for the 
approach to be practical. 

At the present time, there does not exist any analytic non-perturbative theory that is valid for an arbitrary 

roughness. The reason for the lack of such general analytic theory is that for strongly rough surfaces higher order 
scattering processes become important. In consequence the boundary conditions to be satisfied on the random interface 
become dominated by non-local effects. This means that the total field on the surface at some point depends on the 
total field in other locations on the surface. These non-local boundary conditions hamper the development of analytic 
theories for strongly rough surfaces. 

The best one can do at the present for these strongly rough surfaces is to resort to a numerical simulation approach. 
This approach, as we will see below, is based on deriving a set of coupled integral equation for the source functions, 
the field and its normal derivative evaluated on the surface. With the knowledge of these sources, the total field, and 
therefore the solution to the scattering problem, can be obtained from the extinction theorem in any point above the 
surface. We will now outline how all this comes about. 



1. The Extinction Theorem 



We will now derive the so-called Ewald-Oseen extiction theorem first formulated by P. P. Ewald and C. W. Oseen in 
the beginning of this century. The numerical simulation approach to be presented later in this section will be based 
directly upon this theorem. 

Let us start by recall from Sect. IIIB that the primary field, $y(r|a;), satisfies the wave equation 

(dl+E{uj)'^)^,{v\u:) = -JT{v,oj) (111) 

where e{uj) is the dielectric function of the medium where the field is evaluated, 9r = V is the nabla-operator in 
the number of spatial dimensions considered, while J^**(r) is an external source term for the field. In order to solve 

the scattering problem in question, one may solve this equation in the regions of constant dielectric function and 
match these solutions by the boundary conditions that the field, and its normal derivative, should satisfy on any 
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FIG. 7: The geometry considered in the extinction theorem. 



interface (See Sect. II C). However, it is often more convenient to take advantage of certain integral theorems that is 
a consequence of the wave equation (111). The extinction theorem provides such an example. 
The wave equation (111) is accompanied by the following equation for its Green's function [56] 

(d^ + s{uj)'^^G{r\r';uj) = -4nS{r-r'). (112) 

Furthermore, we are only interested in out-going solutions to this equation that fulfill the Sommerfeld's radiation 
condition (at infinity) [107] 

lim r {drG - ikG) = 0, (113) 

r— ^oo 

where r = |r|. 

In two-dimensions, that we consider in the present introduction, an explicit representation of the out-going, free 
space Green's function is provided by [34, 56] 

G(r|r';a;) = (e(,^)^ |r - r'|) , (114) 

where, Hq^\z), is the Hankel-ftmction of the first kind and zeroth-order [39, 56] and r = (xijX^). 

Let us start by considering a spatial region fl containing a homogeneous, isotropic dielectric medium. This region 
has a boundary dSl (See Fig. 7). The exterior of the region fl will be denoted by Q where its boundary is dSl. Notice 
that dtt includes dft in addition to the surface at infinity. We assume that an external source is present somewhere 
in the external region fl and that no sources are present within Q. 

If we multiply Eqs. (Ill) and (112) by respectively G(r|r';a;) and — $i,(r|w), add the resulting equations, and finally 
integrate the result over the exterior region fl we are left with^^ (r' € iT) 

[ dv' [^^{v'\uj)dl.G{v'\v-uj) - dl.^,{v'\uj)G{v'\v-uj)\ 

= -h + { *;(*)• (115) 

Since G(r'|r;a;) is the out-going free space Green's function the first term of the right hand side is just the incident 
field due to the source, i.e. 

^ f dr'jr(r',c.)G(r'|r;a;) = $r(rk). (116) 

This relation holds true independent of r is located in the exterior (fi) or interior (f2) region. 

Furthermore, by taking advantage of Green's second integral identity that for two well-behaved^^ functions u(r) 
and v{r) defined on a region V, reads [34, 56] 

JdV 

(117) 



We have here interchanged r and r' for later eonvenienee. 

By well-behaved we here mean functions that at least are differential two times. 
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where 9„ denotes the outward normal derivative to dV, Eq. (115) can be written as 

$r(rN + ;;^ / dS'[M^'\u;)dn'G{r\v';^)-dn'M^'\Lo)G{v\v'-^)] = j ^ ^ (118) 

where dS' is a surface element. In writing this equation we have explicitly used the fact that the portion of the surface 

integral over dO,, that corresponds to the surface at infinity vanishes due to Sommcrfeld's radiation condition satisfied 
by G(r|r'; ui). Hence the only surface left in the surface integral is dil as indicated in the above equation. In addition 
we have also utilized the relation 9„ = —dn for the outward normal derivative to the region Cl, while dn is the outward 
normal derivative for the same surface, but for region Q. Notice that the incident field term is present due to the face 
that the volume fl contains a source. If this region is source-less this term is missing. 

Eq. (118) with the right-hand-side set to zero is the extinction theorem. It is so named because the incident field is 
extinguished in region Q by the induced field as represented by the second term of the left-hand-side of this equation. 
Furthermore, Eq. (118) with r G fl expresses the fact that the field at any point outside fl can be found by performing 
a surface integral over dfl. In order to do so, however, the total field and its normal derivative on the surface dfl has 
to be known. Hence the scattering problem is equivalent to finding the field and the normal derivative on the surface. 



The Scattered and Transmitted fields 



From the above discussion, we learned that the essential quantities to look for is the field and its normal derivative 
evaluated on the surface. We will now see how these two quantities can be calculated by taking advantage of the 
extinction theorem. This is done by applying Eq. (118) in turn to the different regions naturally defined by the 
scattering geometry as the regions of constant dielectric properties. For the scattering system depict in Fig. 6 this 
means to apply Eq. (118) separately to the regions X3 > Ci(2;i) and X3 < C2(a^i). The result is^^ 

e{xs-axi))^t{r) = €"^(r) + ^ |(ix;7(a;'i)[$^(r')a„'G+(r|r')-5„,$^(r')G+(r|r')]|,,^^(,,), (119a) 

eiCixi) - X3)$;(r) = / dx[ j{x[) [$;(r')5„,G_(r|r') - a„.$;(r')G_(r|r')] |^,^^(^,, . (119b) 

where the superscripts ± indicate solutions to the wave equation (111) in regions of dielectric function s±{ui). Fur- 
thermore, we have defined 

dn = ^^^-^'^-f^^ , (120a) 

where 



7(xi) = \/l + C'(a:i)2. (120b) 

In writing Eqs. (119) we have taken advantage of the assumption made earlier that the surface, C{xi), is a single- valued 
function of xi so that its surface element becomes 

dS = 7(xi)dxi. (121) 

If this assumption does not hold true, the discussion becomes considerably more difficult. A treatment of such a case 
can be found in e.g. Ref. [18]. However, we will not here considered this possibility any further. 

Notice that the integral equations (119) are uncoupled. However, by taking into account the boundary conditions 
to be satisfied on the rough surface X3 = Ci^i), i.e. 

$+(xi,X3;a;)|^^^^(^^) = <i>- {x,,xs:oj)l^^^^^^^^ , (122a) 



dn^tixi,X3; 

(w) 



(122b) 

X3=i(xi) 



In these equations, and some to come, we have suppressed an expUcit reference to the frequency of the incident light in order to make 
the formulae more compact. 
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the two integral equations will be coupled, and Eqs. (119) take on the form 



9{x3-ax,))<^>+ir\Lo) = <i>r(rk) / dx'^ [A+(r|x;;u;)^,KM - B+(r|x'i;u;)AA,KM] , 



6{C{xi)-Xi)^-{v\uj) = - Jdx[ 
where the symbols k,^{oj) have been defined earlier in Eq. (14(i). Here we have introduced the source functions"'^® 



A_{r\x[;u;)Mx[\u;) - ^^^B_{r\x[;w)M4x[\w) 



M^{Xi\u)) = -f{Xi)dn^t{xi,X3\u))\ 



X3=f(xi) 



as well as the kernels 



A±{r\x[;uj) = — jix[) dn'G±{xi,X3\x[,x'ri) 



B±{r\x[;u) = — G±{xi,X3\x[,x'3) 



(124a) 
(124b) 



(125a) 
(125b) 



Notice that the second term on the right-hand-side of Eq. (123a) represents the field scattered from the rough 
surface, $^'^(r|a;). By substituting the following Fourier representation for the Green's function [39] 



G+(r|r';a;) = / ^iq(xi-x',)+ia+(q,w)\x3-x'.,\ 
J_^2'K a+{q,uj) 



(126) 



into Eqs. (125), and the resulting expression into Eq. (123a), we find that the scattered field far above the surface, 
X3 ^ C(^i)i can be written as 

^1\t\uj) = [ ^ i^^(g,a;)e*«^l+^"+(«■'^)^^ (127a) 

where the scattering amplitude is given by the following expression 

RM,^) = TT-f ^ / c/a;i e-^«-^-'"+(«'")«(^^) [i {qC {xi) - a+{q,u)} J^,{xr) - Af^x^)] . 

2a+{q,uj) J_^ 

In these expressions a+{q,Lj), and later to be used a-{q,uj), are defined as in Eqs. (48) and (50). 

If the medium occupying the region .T3 < C(-''^i) transparc^nt, a transmitted field will also exist. It is given by 
the right-hand-side of Eq. (123a). Under this assumption, a Fourier representation for G_(r|r';a;), equivalent the one 
given in Eq. (126), will give a transmitted field in the region X3 <C Ci^i) of the form 



where the transmission amplitude is defined as 
9 1"°° 

TJq,u) = / dxi e-*«^i+*«-(9'"K(^i) 

2a-{q,uj) J_^ 



T^{q,io)e 



iqxi—ia- {q,uj)x3 



(128a) 



i{qC{xi)+a.{q,uj)}Mxi) - ^#\k{xi) 



The Equations for the Source Functions 



Notice that the operator j{xi)dn appearing in J\fv{xi\u!) is nothing else then the unnormalized normal derivative. 
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In order to solve the scattering problem we sec from Eq. (1276) that we need to know the source functions !F^{xi\Ltj) 
and My{xi\ijj). The question therefore is: How to calculate these source functions? A coupled set of equations for 
these sources are most easily obtained by setting X3 — C,{xx) + 77, with 77 — > 0+, in Eqs. (123). Doing so results in the 
following set of inhomogeneous, coupled integral equations for the sources 

T.{xi) = rj^^ixt) Jdx[[A+{x,\x[)J^4x[)-B+{x,\x[)K{x[)], (129a) 





where the kernels are defined as 



= jdx\ 



A-{xi\x[)M<) - ^B-{xi\x[)Af4x[) 



(129b) 



A±ixi\x[) = ^lim A±(r|a;'i)|,^^^(,^)+^, (130a) 
B±ixi\x[) = ^hm S±(r|x;)|,3^^(,3)+^. (130b) 

In order to solve Eqs. (129), the integral equations arc c;onvc!rtcd into matrix equations by discretizing the; spatial 
variables xi and x'^ and using some kind of quadrature scheme for approximating the integrals that they contain. 
First of all, the infinitely long surface is restricted to a finite length Li, so that the spatial integration range from 
— Li/2 to Li/2. Second, a grid defined according to 

= N„ = -^+(''"0^^' n^l,2,3,...,N, (131) 

with — Li/N is introduced for x[. If wc assume that the source functions are slowly varying fmictions over a grid 
ceU (of size A^), they can be considered as constant over this distance and therefore put outside the integral. The 
integral equations (129) are thus converted into the following coupled matrix equations by putting xi = 

N 

MU) = + E [^'innMO - Bt^nJ^^iO] , (132a) 



N r 



(132b) 



where ^^"°(Cm) is defined from Eq. (124a) by using $™'^(a;i, xajw) for the field $+(a;i, xalw). Moreover, the matrix 
elements ^"^^ ^mn defined as 

= / dx[A±{UW,), (133a) 

■/€„-A€/2 
/.|„+A5/2 

B±„ = / dx[B±iUx[). (133b) 

■/€„-A€/2 

It should be kept in mind that these matrix elements are related to the Hankel function, Hq^^ (z), and its derivative, 
through the (two-dimensional) Green's hmction that enters via Eqs. (125) and (130). Care has to be taken when 
evaluating these matrix elements since the Hankel functions are singular when their arguments vanish. Hence the 
kernels, ^±(a;i|a;'i) and S±(xi|x'i), are also singular when xi = x[. Fortunately these singularities are integrable 
so that the matrix elements, and in contrast to the kernels, are well define everywhere. The somewhat 

technical procedure for showing this is presented in Appendix A, from where we obtain that (see Eqs. (^11) and 
(A12)) 

, f A^A±(C„|^„), m^n, 

Amn = i 1 , ("M (134a) 



2 + '^^ 4^72(4. 



m 



and 
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The matrix equations (132), together with the expressions for the matrix elements Eqs. (134). can readily be 
put onto the computer and solved by standard techniques from linear algebra [40, 41] in order to obtain the source 
functions. With these source functions available, the scattering amplitude, and, if defined, the transmission amplitude, 
can be obtained from respectively Eqs. (1276) and (1286). These amplitudes are again related to physical observable 
quantities, like the mean differential reflection or transmission coefficients, as discussed earlier. Hence the scattering 
problem is in principle solved! 

It should also be mentioned that the approach presented here can be generalized to more complicated scattering 
geometries like film systems etc. [22, 68, 122]. However, in such cases the higher demand is put on computational 
resources. 

2. A Remark on the Accuracy of the Numerical Simulation Approach 

The numerical approach described above is formally exact since no approximations have been introduced. It is 
therefore in principle applicable to scattering from surfaces of any roughness. It has proven useful in many situations, 
and serve today as a standard, and invaluable, tool for rough surface scattering stiidies. This is in particular true for 
scattering from strongly rough surfaces where it represents the only rigorous available method at present. 

Even if this approach is formally exact, it has some practical limitations. Imagine a weakly rough metal surface 
that is illuminated by p-polarized light. In this case the incident light can excite surface plasmon polaritons that will 
travel along the rough surface. The mean free path of these surface plasmon polaritons will in the present case be 
quite large. 

In computer simulations we are, of course, not able to represent infinitely long surfaces. Instead we are limited to 
surface of finite length. To avoid essential contributions to the simulation results from artificial scattering processes 
e.g. where surface plasmon polaritons are being scattered from the edges of our (finite length) surface, its length 
needs to be long. In order not to compromise the spatial resolution used in the simulations, big demands on computer 
memory and CPU-time is a consequence. This sets a practical limit for the \isc of rigorous numerical simulations 
for weakly rough surfaces. However, for such kind of roughness, perturbation theory, where we by construction are 
using surfaces of infinite length, arc adequate and accurate as discussed earlier. The present limitation of the rigorous 
numerical simulation approach should therefore not represent a too severe restriction from a practical point of view. 

We should also mention that there exists another numerical technique that for the same amount of memory used 
in the rigorous approach can handle much longer surfaces (and therefore reduce edge effects) . This technique is based 
on a numerical solution of the reduced Rayleigh equation, that was introduced in Eq. (78) of Sect. Ill E, and is the 
single integral equation satisfied by the scattering amplitude. The interested reader should consult Refs. [42] and 
[121] for details regarding this numerical technique. This approach is obviously restricted to surfaces for which the 
Rayleigh hypothesis is valid. It can therefore not be applied to surfaces of arbitrary roughness, but it is valid for 
surfaces that practically can not be treated within perturbation theory. A direct numerical solution of the reduced 
Rayleigh (integral) equation can therefore be looked upon as a bridge between perturbation theory and the rigorous 
numerical simulation approach. 

IV. PHYSICAL PHENOMENA IN ELECTROMAGNETIC ROUGH SURFACE SCATTERING 

Wave scattering from randomly rough surfaces has a long history in science [1-4]. In the overall majority of 
theoretical studies conducted up to the early 1980 's, single-scattering approaches were used [5-8]. However, around 
this time people started getting interested in the effects and consequences of incorporating multiple-scattering events 
into the theories. It created a lot of excitement in the field when new and interesting multiple scattering phenomena 
were either predicted theoretically and/or observed in experiments. During the period of time that has passed since the 
early 1980's, multiple scattering effects from randomly rough surfaces have attracted much attention from theorists 
and experimentalists alike, and today the research in this field is concentrated around different kinds of multiple 
scattering effects [22]. 

It ought to be mentioned that many of the effects to be discussed here are not exclusive to surface scattering. Quite 
a few of them have in fact their analogies in light scattering from volume disordered systems. For discussion of light 
scattering from volume disordered systems the interested reader is referred to the literature [59, 60]. 

In this section we aim at discussing some of the new multiple scattering effects that might take place when electro- 
magnetic waves are scattered from a randomly rough surface. The technical details on which the present section rely 
were mainly presented in the previous section. We have therefore tried to keep the discussion at a phenomenological, 
and hopefully pedagogical, level. Unnecessary technical details have been avoided whenever possible. As a service to 
the more technical oriented reader, an extensive reference to the original literature has been made. 
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FIG. 8: Pcrturbativc calculations for the mean difFerential reflection coefficient for the incoherent component of the light 
scattered from a randomly rough silver surface. The incident angles of the light of wavelength A = 632. 8nm were (a) ^0 = 0° 
and (b) On = 25°. The dielectric constant of silver at this wavelength is e{ijj) = —7.5 + iO.24. The surface was characterized by 
a Gaussian height distribution of RMS-height 5 = 5nm and a Gaussian height-height correlation function of correlation length 
a = lOOnm. 



A. Coherent Effects in Multiple-Scattered Fields: Weak Localization of Light on a Randomly Rough Surface 

In 1985 McGurn, Maradudin, and Colli [54] predicted theoretically the existence of what later has been known as 
the enhanced backscattering phenomenon in surface scattering. This phenomenon express itself as a well-defined peak 
in the retroreflection direction of the angular dependence of the light scattered incoherently from a rough surface. 
The work in Ref. [54] was the first to report on an effect that was shown to be caused by multiple scattering processes 
taking place at the rough surface. The enhanced backscattering phenomenon is an example of what is known as a 
coherent effect in the multiple scattered field. Later on, other coherent phenomena, like the enhanced transmission [61], 
satellite peaks [65, 66] and enhancements due to the excitations of magnetoplasmons [67, 68] were predicted. 

1. Enhanced backscattering 

In this subsection the backscattering enhancement phenomenon will be discussed. Since the mechanisms that give 
rise to it are different for weakly and strongly rough surfaces, they will be treated separately. The scattering system 
that will be considered is depicted in Fig. G and consists of a single rough vacuum-metal surface. 

Weakly Rough Surfaces 

It is familiar from every day life that if the surface is not too rough, the waves incident on it will mostly be scattered 
into the specular direction. That is, if the angle of incidence is 9o, then most of the energy will be scattered into the 
direction 6s = Oo, which defines the specular direction. For a weakly rough surfaces the intensity, or equivalently the 
mean differential reflection coefficient (DRC), will have a maximum a specular peak for the scattering angle 
= ^0- Normally the specular peak is not of any interest to us, and it is therefore in theoretical studies usually 
subtracted of, leaving only the intensity that results from light scattered incoherently by the rough surface. 

In 1985 McGurn, Maradudin, and Celli [54] predicted based on a perturbation theoretical study, that also in the 
retroreflection (anti-specular) direction of the angular dependence of the mean DRC there might be an enhancement. 
This effect, known today as enhanced backscattering. manifest itself as a well-defined peak in the retrorefiection 
direction of the angular dependence of the intensity of the light that has been scattered incoherently from the random 
surface. 

In the original paper by McGurn et al. [54] , the calculation of the enhanced backscattering peak was carried out for 
p-polarized light scattered from a weakly rough silver surface. Their calculations, based on a many-body perturbation 
theory, took into account multiple scattering events in the calculation of the intensity scattered incoherently by the 
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FIG. 9: Diagrams showing two of the scattering events that through interference give rise to the enhanced backscattering peak 
for weakly rough surfaces. 



surface. In Figs. 8 we show the results of a small amplitude perturbative calculation, like the one given in Sect. IIIF, 
for the incoherent contribution to the mean differential reflection coefficient for p-polarized light incident at angles 
^0 = 0° (Fig- 8a) and = 25° (Fig. 8b) on the rough silver surface. Terms to 4'th order in the surface profile 
function were included which is enough to include all double scattering processes. The wavelength of the incident 
light was A = 632. 8nm, and the dielectric constant of silver at this wavelength is e{u) = — 7.5 + i0.24. The surface was 
assumed to be characterized by a Gaussian height distribution and the height-height correlation function was also of 
the Gaussian type. The root-mean-square (rms) height of the surface was d = 5nm while the correlation length was 
a = lOOnm. Prom Fig. 8 we see well-pronounced peaks for the retroreflection directions. It should be stressed that 
it is the incoherent component of the mean DRC that is plotted in these figures, so that the peak seen at Og = 0° in 
Fig. 8a is no specular effect^" since all contributions from specular scattering have been subtracted off. That this is 
the case should be obvious from the position of the enhanced backscattering peak seen in Fig. 8b corresponding to 
the incident angle 9o = 25°. 

The natural question is now: What is the origin of the enhanced backscattering peak? It was realized that it had 
to be caused by multiple scattering since it had not been seen earlier when using single scattering theories. It turned 
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FIG. 10: The same as Fig. 8a, but now showing separately the 2nd and 4th order contribution in the surface profile function 
to the mean differential reflection cocflicicnt. The sum of those terms gives the curve in Fig. 8a. Notice that the enhanced 
backscattering peak comes from the 4th order contribution, i.e. form the double scattering contribution. 
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Recall that for normal incidence the specular and anti-specular directions coincide. 
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FIG. 11: Experimental results (open circles) for {dRp/dOs} ^^^^i^ as a function of the scattering angle 0s for three different 
incident angles when p-polarizcd light of wavelength A = 612. 7nm is incident on a onc-dimcnsional random gold surface 
{e{uj) = —9.00 + il.29). For the random surface a West-O'Donnell power spectrum with fc_ — 0.82(a;/c) and k+ = 1.29(a'/c) 
was used. The RMS-height of the surface was 5 = 10.9nm. The solid lines and the open triangles are perturbation theoretical 
results based on respectively small amplitude and many-body perturbation theory. (After Ref. [20]). 



out that the origin lies in the interference between a multiple scattered path with its reciprocal partner [20, 54]. To 
illustrate this, let us consider the double scattering path shown in Fig. 9a. Here an incident wave excites through the 
breakdown of infinitesimal translation invariance of the system, a surface plasmon polariton that propagates along 
the surface. At the next scattering event this surface polariton is converted back into a volume electromagnetic wave 
that propagating away from the surface. This path has a reciprocal partner (Fig. 9b) where the scattering takes 
place from the same scattering centers at the rough surface, but now in the opposite order. For the backscattering 
direction these two paths will have exactly the same amplitude and phase, i.e. they will be coherent, and hence they 
will interfere constructively. However, as we move away from the backscattering direction, the two paths fast become 
incoherent so that their intensities just add. Thus, due to the interference nature of the enhanced backscattering peak 
the amplitudes at the position of the peak would in the absence of single scattering be twice that of its background due 
to the cross-terms originating from the square modulus of the amplitudes needed in order to calculate the intensity. 
However, notice that it is not uncommon that single scattering gives considerable contribution to the mean differential 
reflection coefflcient of the light scattered incoherently from the surface. In such cases, the height of the peak is not 
twice of its background. 

To show that multiple scattering indeed is the origin of the enhanced backscattering phenomenon, we show in 
Fig. 10 the different contributions to the incoherent component of the mean DRC obtained from Eq. (89a). We recall 
that the first term of this equation is the single scattering contribution, i.e. it is of 2nd order in the surface profile 
function Ci^i)- The next two terms are both double scattering contributions or, cquivalcntly, 4th order contributions 
in the surface profile function. From Fig. 10 it is seen that the single scattering contribution (2nd order in Ci^i)) is 
a smooth function of the scattering angle. Furthermore, it is seen that the peak stems from the double scattering 
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FIG. 12: Numerical Monte Carlo simulation results for the mean difFcrcntial reflection coeflicient, (,dRv /dOs) ^^^^f^, for (a) p- 
and (b) s-polarized light scattered from a rough silver surface of RMS-liciglit 5 = 15nm. The angle of incidence was ^0 = 0° 
and the wavelength of the incident light was A = 632. 8nm. At this wavelength the dielectric constant of silver is e(u)) = 
—7.5 + iO.24. Furthermore, the surface was of the West-O'Donnell type characterized by the parameters k- = 0.82(a;/c) and 
k+ = 1.92(a;/c). The simulations were performed by numerically solving the reduced Rayleigh equation that the scattering 
amplitude satisfies [42, 121]. The number of samples used in obtain these results was A^^ = 3000. Notice that in the case of 
s-polarization there is no enhanced backscattering peak in contrast to what is the case for p-polarization. This difference is 
caused by the fact that s-polarized incident light can not excite surface plasmon polaritons at a rough vacuum-metal interface. 



contribution, i.e. it comes from the second and third terms of Eq. (89a). In a diagrammatic language this term comes 
from the maximally crossed diagrams. The interested reader should consult Ref. [54] for additional details. 

It should be noted that even if we earlier only included fully the lowest order multiple scattering process (double 
scattering), higher order processes will not change the statement that the enhanced backscattering phenomenon is 
caused by multiple scattering through the constructive interference between a scattering path with its reciprocal 
partner. 

The enhanced backscattering effect from weakly rough vacuum-metal surfaces was observed in experiments by West 
and O'Donnell [20] in 1995 in the scattering of p-polarization light from a rough gold surface of RMS-height 5 = 10.9nm. 
The power-spectrum used in these experiments was of the rectangular type also known as the Wcst-O'Donnell power- 
spectrum. The remaining parameters used are defined in the caption of Figs. 11. We have in Figs. 11 reproduced 
their experimental results (open circles) together with some perturbation theoretically results (solid lines and open 
triangles). At least for the two smallest angles of incidence well-defined peaks around the retroreflection direction in 
the experimental results are seen. 

For weakly rough surfaces we just argued that the origin of the enhanced backscattering effect involves surface 
plasmon polaritons. In s-polarization, a rough (one-dimensional) vacuum-metal interface does not support such 
surface waves. Hence, one does not expect to see any backscattering peak for this polarization for weakly rough 
surfaces. This is indeed seen from Figs. 12 showing numerical simulation results for p- and s-polarized incident light 
based on the solution of the reduced Rayleigh equation that the scattering amplitude satisfies. The power spectrum 
used for the surface was again of the West-O'Donnell type, and it was defined by the parameters fc_ = 0.82(a;/c) 
and /c+ = 1.92(a;/c). It is seen from Figs. 12 that only in p-polarization do we sec an enhanced backscattering peak. 
Prom Eq. (89) we see that the single scattering contribution is proportional to the power spectrum, (7(|A;|), of the 
surface roughness. Hence, if fc_ > the incoherent component to the mean DRC, {dRu/dOs)^^^^}^ should not contain 
any contribution from single scattering events in the angular range —9- < 6s < 0-, where 0- = arcsin ((fc_c)/a;). 
With the parameters used in obtaining Figs. 12 this gives 0- = 55.1°. For scattering angels 16*5 1 > 0- = 55.1° single 
scattering is allowed. This can be seen as a jump in Figs. 12 at this angle. Furthermore, around the backscattering 
peak, single scattering should be absent and indeed the enhanced backscattering peak is twice that of its background 
as predicted above. Notice also that the overall fraction of the light scatterer incoherently from the surface is at least 
one order of magnitude lower for s- then p-polarization. 



Strongly rough surfaces 



43 




-30 30 
93 (deg) 



60 90 





-30 30 
03 (deg) 



60 90 



-30 30 60 90 
63 (deg) 




-30 30 60 90 
63 (deg) 



FIG. 13: Rigorous Monte Carlo simulation results for the angular dependence of the incoherent coniponont of the mean 
differential reflection coefficient for the scattering of (a and b) p- and (c and d) s-polarized incident light from a strongly 
rough silver surface. The wavelength of the incident light was A = 612. 7nm for which the dielectric constant of silver is 
e(w) = —17.2 + iO.498. The incident angles of the light were (a and c) ^0 = 0° and (b and d) 60 = 25°. The strongly rough 
surface was characterized by Gaussian height distribution of RMS-height 5 = 1.2/xm and the transverse correlation length for 
the Gaussian correlation surface was a = 2/Ltm. The length of the surface was L = 25.6/xm and a finite sized beam of width 
g = 6.4/im was used in the simulations. The number of discretization points was = 500. The numerical results were all 
based on an ensemble average over = 3000 realizations of the randomly rough surface. 



We will now consider strongly rough surfaces. In order to study the backscattering phenomenon for such surfaces we 
have to resort to numerical simulations like e.g. the approach outlined in Sect. Ill I. In Figs. 13 we present the results 
of such simulations for the angular dependence of the incoherent component of the mean DRC for (a and b) p- and (c 
and d) s-polarizcd light incident on a rough vacuum-metal surface of RMS-hcight S = 1.2/im. The correlation length 
for the Gaussian correlated surface was a = 2/im. The main difference between these results and those for the weakly 
rough surfaces presented earlier (Fig. 12) is that we now also observe an enhanced backscattering peak in the case of 
s-polarization. So what is the reason for this difference between weakly and strongly rough surfaces when it comes to 
the backscattering phenomenon? The explanation lies in the mechanism causing the backscattering peak for strongly 
rough surfaces [21, 62-64]. Since the excitation of surface plasmon polaritons is weak for strongly rough surface, it is 
unlikely that the reason for the backscattering peak is caused by this type of surface waves. Such a mechanism could 
not in any case explain the presence of the backscattering peak observed for s-polarization. Instead the backscattering 
peak for strongly rough surfaces arises due to the constructive interference between multiple scattered volume paths 
like e.g. those shown in Figs. 14. In this case no surface waves are excited, but instead the multiple scattering takes 
place within the valleys of the now strongly rough surface. The incident wave, that after its first encounter with 
the rough surface, is scattered at least one more time before leaving the surface for good. Also in this case for the 
backscattering direction this path has a reciprocal partner that is phase coherent with the first one and with which 
the latter path can interfere constructively. Since this mechanism does not involve any surface plasmon polaritons, 
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FIG. 14: Diagrams showing two double scattering paths that for strongly rough surlaces through interference represent the 
main contributes to the enhanced backscattering peak phenomenon. 



there is no reason why the backscattering phenomenon should not show up also in s-polarization from strongly rough 
surfaces. In fact as can be seen from Fig. 13c and d, the backscattering peak in s-polarization is as pronounced as for 
p-polarization. Observe also that the energy scattered incoherently, which for strongly rough surfaces is close to the 
total scattered energy, is of roughly the same order for both polarizations. This is in contrast to the situation found 
for weakly rough surfaces. 

The enhanced backscattering phenomenon from strongly rough surfaces was experimental confirmed as early as 
1987 by Mendez and O'Donnell [62]. This was just two years after its theoretically prediction by McGurn et al. [54] 
for weakly rough surfaces. In fact these experiments provided the first experimental evidence what-so-ever for the 
enhanced backscattering phenomenon from rough surfaces. 



2. Satellite Peaks 



The backscattering phenomenon discussed in the previous subsection is not the only coherent effect that might 
exist when light is scattered from a randomly rough surface. Another such effect is the existence of so-called satellite 
peaks predicted by Freilikher, Pustilnik, and Yurkevich [65] in 1994. Satellite peaks are enhancements in the angular 
distribution of light scattered incoherently from scattering systems that supports more the one surface [65-68] or guided 
waves [42, 69-71, 109, 121]. As will be shown in detail in Subsect. IV A3, they are not caused by interference between 
reciprocal paths as was the case for the backscattering phenomenon, but instead by interference of nonreciprocal paths. 
These enhancements should occur for scattering angles that are located symmetrically with respect to the position of 
the enhanced backscattering peaks that the scattering system also gives rise to. 

To illustrate this, let us study the film scattering system shown in Fig. 15. Here the lower interface is rough and 
the upper one is planar. Furthermore, the lower semi-infinite medium is assumed to be a perfect conductor, while 
the incident medium is assumed to be vacuum. In Figs. 16 we show the results of numerical simulations for the mean 
differential reflection coefficient in the case of s- (Fig. 16a) and p-polarizcd (Fig. 16b) incident light. The remaining 
surface parameters arc given in the caption of Figs. 16. In the case of s-polarization, the scattering system of mean 
thickness d = 500nm supports two satellite peaks, while in the case of p-polarization it can at most support six such 
peaks [121]. The positions of these peaks are indicated by dashed vertical lines in Figs. 16. From Fig. 16a the two 
satellite peaks that the scattering system supports in this case are easily distinguished from the background. However, 
from Fig. 16b one sees that only four out of the six possible peaks can be observed. There are two reasons why some 
of these satellite peaks may not be observable: First, some of them may lie in the non-radiative part of the spectrum, 
and are therefore not even in principle observable. Second, their strength might be to low to be observable [42, 121], 
i.e. one (or both) of the channels involved in the interference process that gives rise to the satellite peaks might have 
too low intensity [42, 121] (see Subsect. IV A 3). 

It can be shown (result not shown) that by reducing the thickness of the film, and thus reducing the number of 
guided waves that the system supports say to one, all the satellite peaks vanish while the enhanced backscattering 
peak is still present [71]. In an analogous way, if the film thickness is increased, more then two satellite peaks might 
be seen [71]. 



45 




FIG. 15: A sketch of a film scattering geometry that supports guided waves and that may give rise to satellite peaks in the 
angular dependence of the scattered light. 



3. A Formal Approach to Enhanced Backscattering and Satellite Peaks 



In the previous two subsections the enhanced backscattering and satellite peaks phenomena were discussed. In the 
present subsection a more detailed analysis and formal approach towards these two phenomena will be presented. In 
particular we will determine at which positions the satellite peaks arc to be expected. 

Let us consider a general film scattering system, where at least one of the interfaces are rough. Fig. 15 provides one 
example of such a system. Depending on the thickness of the film, the scattering system supports > guided waves 
at the frequency u of the incident light. The wavenumbers of these modes, or "channels" as some authors prefer to 




FIG. 16: The contribution to the mean differential reflection coefficient from the incoherent component of the scattered fight 
{dRv / d9)^^^^j^ as a function of the scattering angle 6s when an s- (Fig. 16a) or p-polarizcd (Fig. 16b) plane wave of wavelength 
A = 633nm is incident normally (Oq =0°) on the film scattering geometry shown in Fig. 15. The dielectric constant of the film 
at the wavelength of the incident light is e{uj) = 2.6896 + iO. 01, and the films mean thickness was d = 500nm. The semi-infinite 
medium which the film is ruled on was a perfect conductor. The surface profile function C,{x\) of the film-conductor interface 
was characterized by a Gaussian surface height distribution of RMS-height 5 = 30nm and a West-O'Donnell power spectrum 
defined by the parameters fe_ = 0.82(aj/c) and k+ — 1.97(a;/c). The length of the surface used in the simulations was L = 160A. 
The dashed vertical linos indicates the estimated positions of the satellite peaks (sec Ref. [121] for details). The results were 
obtained by numerical simulations based on the reduced Rayleigh equation. The data in Fig. 16b have bee smoothed to make 
the positions of the satellite peaks more apparent. (After Ref. [121].) 
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FIG. 17: Illustration of two double scattering sequences occurring in the scattering of electromagnetic waves from a bounded 
systems that supports more the one guided (or surface) wave. 



call them, will be denote by where n = 1, . . . N. Through the surface roughness the incident light may couple 

to these guided waves. 

In Fig. 17 we show two general double scattering paths^^ where the scattering takes place at the same scattering 
centers, but in the reverse order. Such paths will in general be phase incoherent due to the randomness of the rough 
surface. However, wc will now looking into if there arc particular angles of incidence and scattering for which these 
two paths are phase coherent. Let us start by assume that path ABCD goes through channel m and path ACBD 
through the n-channel. The phase difference between these two paths can then be expressed as 

^(t>nm = rsc ■ (ko + K) + \rBc\ - qrn{^)] ■ (135) 

Here ko and kg are the wave vectors of the incident and scattered waves, respectively, while tbc is the distance 
(vector) from point B to C. According to its definition, we will have phase coherence when this phase- difference is 
zero, i.e. when A(f>nm = 0. Now let us consider separately two cases: (i) n = m and (ii) n ^ m. In the first case 
the last term in Eq. (135) is zero with the consequence that one has phase coherence if kg = —kg. This coherence is 
obviously what gives rise to the enhanced backscattering phenomenon. In the second case when m ^ n, the last part 
of Eq. (135) does not vanish. The condition for phase coherence then becomes 

1 c 

sinOs = -sin6>o± , - |g»('^) - gm(t^)| • (136) 
Veo('^) ^ 

In this equation we have also allowed for the case m = n since it naturally includes the position of the backscattering 
peaks. Hence, Eq. ( 136) defines the angles for which peaks due to coherent effects arc expected in the angular 
dependence of the light scattered incoherently from the randomly rough surface. The angle obtained for m = n is the 
position of the backscattering peak, while the angles obtained for m ^ n correspond to satellite peaks. The reader 
should check that the angles obtained from Eq. (136) fit the position of the satellite peaks shown in Fig. 16. The 
values for qn{oj) can be found in Rcfs. [42] and [121]. 

This concludes our discussion of coherent effects in the scattered field. Even though we have focused on the reflected 
light, it should be pointed out that there also exist similar effects in the transmitted light [61]. For a discussion of 
this case the reader is referred to the literature for details [22, 61]. 

B. Localization 

The notion of localization was introduced into physics by P.W. Anderson in his famous 1958 paper [72], a work 

that he was awarded the Nobel Prize for. Anderson studied the transport properties of electrons in materials with 
bulk disorder. This study led him to what today is known as the Anderson localization phenomenon, sometimes also 



Wc here consider double scattering for simplicity. Higher order scattering processes can be treated the same way, but doing so will not 
change the main conclusions. 



47 




called strong localization. The phenomenon expresses itself by a disorder-induced phase transition in the transport 
behavior of the electrons. As the disorder is increased in a three dimensional system, the scattering evolves from a 
diffusion regime, for which the well-known Ohm's law holds, to a localized regime in which the material behaves as 
an insulator and all states are localized in space. These two phases are separated by the mobility edge. Anderson 
suggested [72] that the localization regime was caused by strong interference that resulted in an exponential decay 
of the wave function of the electrons in all directions and a subsequent vanishment of the diffusion constant. Hence 
localization is a miiltiple scattering phenomenon. In contrast to what is the case for three dimensional systems, all 
states are expected to be localized for systems that are one- and two-dimensional [73] . However, in this latter case it 
might happen that the localization length is large, and even exceeds the sample size. For a more detailed introduction 
to localization the reader is directed to Rcfs. [73] and [74]. 

It was realized shortly after Anderson published his pioneering work [72], that a similar phenomenon should also 
be expected for multiple-scattering of electromagnetic waves. At room temperature, photons can be treaded as non- 
interacting. They are therefore not hampered by the troublesome self-interaction that electrons have and that are 
known to represent another, but different mechanism towards an insulator regime [75]. This fact makes photon dis- 
ordered systems ideal for studying Anderson localization [76]. However, it should still take several decades before 
localization of electromagnetic waves was confirmed experimentally. Finally in 1997 Wiersma, Bartolini, Lagendijk, 
and Righini were able to obtain direct experimental evidence that confirmed the localization hypothesis for electro- 
magnetic waves in disordered media [77]. These experiments were performed on a system containing very strongly 
scattering semiconductor powders. Thus, the scattering system involved was of the bulk disordered type. 

It is still, however, an open question if Anderson localization of electromagnetic waves can be observed experimen- 
tally for systems containing only surface disorder, even though it should exist in principle due to the system being 
two-dimensional (and in some cases effectively one-dimensional) . Since there is only disorder on the surface, localiza- 
tion can only exist for electromagnetic waves that happen to "live" on or close to the surface. Such waves are called 
surface waves, and we will here focus on surface plasmon polaritons (SPP) that might exist on e.g. a vacuum-metal 
interface. SPP localization should be characterized by the exponential decay of the transmitted intensity as a function 
of distance traveled by the SPP along the rough surface. 

However, the problem of observing SPP localization for surface disordered systems is that Anderson localization 
might be masked by more dominating effects giving rise to the same type of signature as localization itself — the decay 
exp {—L/£t{<-^)) with system size L of the transmittancc where ixi^) is a decay length. The competing effects arc in 
addition to the Anderson localization: (i) ohmic losses in the metal due to Ime(uj) ^ 0, and (n) roughness-induced 
conversion of surface plasmon polaritons into volume waves above the surface — so-called leakage. Hence the decay 
length, It{'-^)i of the transmission coefficient (that we will define below), should be related to the decay length due 
to ohmic losses, (.^i'^), the one due to leakage, iradi^^), ^"^^ t;he Anderson locahzation length £{uj), according to the 
formulae 



In order to determine the Anderson localization length one has to sort out the contribution from each of these 
competing effects. In other words, we have to identify the lengths £e(a;) and i!rod('^) in order to be able to estimate 

The decay rate due to ohmic losses is easily determined and doesn't represent any serious problem (see discussion 
below). However, a much more severe problem is how to separate the contribution from leakage and localization. 
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Leakage is expected to be a rather strong effect with the consequence that £rad{^) is smah compared to the other 
lengths appearing on the right-hand-side of Eq. (137). If this is the case, a measurement of the transmission coefficient 
for the SPP will result in irii^) — ^rad(t^)- It is therefore important, if we are trying to estimate £{uj), to be able to 
separate the localization length from the one of leakage, or to be able to suppress leakage. The approach we will follow 
here is the latter one — the suppression of leakage. This can be done by specially designing surfaces that suppress 
leakage. How this can be done will be presented briefly below (see Ref. [118] for more details). 



1. The Scattering System 

The scattering system that will be considered in this section is depicted in Fig. 18. We study the scattering of a 

p-polarized surface plasmon polariton of frequency lu propagating in the positive xi-direction and is incident onto a 
segment of a one-dimensional randomly rough surface defined by the equation xs = C{xi)- The surface profile function 
(^{x\) is assumed to be a single-valued function of Xi that is nonzero only in the interval —L/2 < xi < L/2. The region 
.T3 > C,{xi) is vacuum; the region 2:3 < C(-'^i) is a metal characterized by an isotropic, frequency-dependent, complex 
dielectric function s{ijj) = £i{uj) -\- ie2{oj)- We are interested in the frequency range in which £i{uj) < — l,£2(w) > 0, 
within which surface plasmon polaritons exist. Furthermore, the rough portion of the surface is assumed to constitute 
a Gaussian random process and with the other "standard" properties described in Sect. II F. 



2. Surfaces that Suppress Leakage 

The first step towards the estimation of the Anderson localization length for this scattering system is to construct 
surfaces that suppress leakage. We recall that the incident surface plasmon polariton has a wave vector given by 

= fci(w) +ifc2(w), (138) 

where ki (w) and ^2 (u)) are the real and imaginary part of the (complex) SPP wave vector and deflned explicitly in 
Ref. [118]. 

By interaction with the surface roughness, the incident SPP picks up momenta available in the power-spectrum of 

the roughness, and due to scattering, changes its wave vector into q. If this momenta satisfies \q\ < oj/c, leakage has 
occurred. To prevent this, or at least reduce this effect, we might use an intelligently choice for the power-spectrum. 
Such a choice is provided by a (rectangular) West-O'Donnell power-spectrum of width 2Ak located around ±fci(a;) (see 
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FIG. 19: A sketch of the power-spectrum used in order to suppress leakage. See text for details 
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FIG. 20: Numerical calculations for (yijc)^ <^|i?(g|fc)|^^ as a function of cg/tj for a rough silver surface characterized by the 
parameters Afc = 0.3(fci(tj) — (ijj/c)) and 5 = 30nm. The rough portion of the surface had length L — 20A. The frequency 
of the surface plasmon polariton, k{ui) — ki{uj) + ik2{u;) = (1.0741 + i0.0026)tj/c, corresponds to a vacuum wavelength of 
A = 457. 9nm, and the dielectric constant of silver at this frequency was e{uj) = —7.5 + iO.24. The results for 50 realizations of 
the surface profile function were averaged to obtain the results plotted in this figure. 



Fig. 19), i.e. a power-spectrum of the form 

g{\Q\) = - k-)e{k+ - g) + 6{-q - k_)e{k+ + q)], (139) 

where k± = 2ki{uj) ± Ak and Afc must satisfy the inequality Ak < ki{Lu) — — . 

That a surface characterized by the power spectrum (139) suppresses leakage can be seen from the following argu- 
ment: The incident surface plasmon polariton has a wave number whose real part is ki{uj). After its first interaction 
with the surface roughness the real part of its wave number will lie in the two intervals (3fci(a;) — Ak, 3ki{u!) + Ak) 
and (— fci(a;) — Afc, —ki{uj) + Ak). For the same reason, after its second interaction with the surface roughness the real 
part of the wave number of the surface plasmon polariton will lie in the three intervals (5fci(a;) — 2 Afc, 5fci(w) 4- 2 Afc), 
(fci(aj) — 2 Afc, fci(a;) -f 2Afc), and (— 3fci(w) — 2 Afc, — 3fci(a;) -|-2Afc). After three interactions with the surface roughness 
the real part of its wave number will lie in the four intervals (7fci(a;) — 3Afc, 7fci(w)-|-3Afc), (3fci(w)— 3Afc, 3fci(u;)-|-3Afc), 
(— fci(ijj) — 3Afc, —ki{Lu)+3Ak), and (— 5fci(ijj) — 3Afc, — 5fci(ti') + 3Afc), and so on. Thus, for example, if — fci(aj)-|-3Afc < 
— (w/c), so that Afc < i(fci(ti;) — (lu/c)), after three scattering processes the surface plasmon polariton will not have 
been converted into volume electromagnetic waves. In general, if we wish the surface plasmon polariton to be scattered 
n times from the surface roughness without being converted into volume electromagnetic waves, we must require that 

Ak<-(ki(uj)-(uj/c)). (140) 
n 

To illustrate that the above procedure really works, we present in Fig. 20 numerical simulation results for the 
scattering ampHtude above the surface, [iojcf' (|i?^((7, w)P) as a function of cqjio for a silver surface characterized by 
the parameters Afc = 0.3(fci(w) — (w/c)) and 5 = 30nm. This surface should thus suppress leakage up to and including 
third order scattering processes. We see from Fig. 20 that (|i?^(5,tj)p) indeed becomes heavily suppressed in the 
range — (w/c) < q < {(^/c). Notice that the six peaks seen in Fig. 20 correspond to the real parts of the wavenumbers 
of the scattered surface plasmon polaritons resulting from the scattering of an incident surface plasmon polariton of 
wave vector k{uj) = fci(cj) -t- ik2{uj) = (1.0741 -t- z0.0026)a;/c. 
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FIG. 21: Numerical simulation results for (ln|l + t(tj)j^^ vs. the length ol the rough portion ol the surface L. The remaining 
surface parameters were the same used in Fig. 20. The number of realizations used for the rough surface in calculating the 
ensemble average was for each system size of the order of 10^. The error-bars indicate the spread in |1 + t{uj)f that was part 
of this average. The solid line, which slope is related to the Anderson localization length, is a x^-fit to the numerical data 

corresponding to a localization length of £{ui) = (5319 ± 905) A. The length due to ohmic losses in the metal is ^^(tj) = 30.1 A. 



3. The Anderson Localization Length for Surface Plasmon Polaritons Localized on a Randomly Rough Surface 



The transmission coefficient for surface plasmon polariton, T(L), is defined as the fraction of the energy flux entering 
the random segment of the metal surface at xi — —L/2 and that leaves it at xi = L/2, i.e. 



T{L) = 



^inc (^-j) 



where Pinc(2;i) and Ptr(a;i) denote the incident and transmitted flux at position xi. 
Above the surface the field can be written as [118] 



H'^{xi,X'i\ui) 



dq 
2^ 



t{uj)e 



, xi > L/2, 



with Poito) — iaQ{ksp{uj),LL!) and 



£2(tj) - 1 

Thus the transmission coefficient can alternatively be written as [118 

T[L) = |l+t(tj)|'exp 



L 



(141) 



(142a) 



(142b) 



(143) 



where (e{uj) — l/{2k2{u!)) is the SPP mean free path due to ohmic losses. Notice that Eq. (143) separates the 
contribution due to ohmic losses from the one of Anderson localization (and leakage). 
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The quantity that we will be interested in when studying the localization phenomenon is not the transmission 
coefficient itself, but instead hiT[L) or more precisely its (ensemble) average given by 

(InT(L)) = (ln|l+i(a;)|^)-^. (144) 

If the effect of leakage can be neglected, the first term on the right- hand-side of the above equation will only incorporate 
contributions from (Anderson) localization. One therefore writes 

(ln\l + t{io)f) = const.- (145) 

Hence, under the assumption that leakage can be neglected the localization length, ^(w), can be obtained from a 
straight line fit to ^ln|l + f(a;)|^^ vs. system size L. 

In Fig. 21 we present such a plot resulting from numerical simulations using the same surface parameters that lead 
to the results shown in Fig. 20 except that now the length of the rough portion of the surface was different. These 
simulations were performed on the basis of the reduced Rayleigh equation that T{q,w) satisfies [118], 

T{p,u) = Vip\ksp{uj)) + J ^V{p\q)Go{q,u;)T{q,u), (146) 

where V{p\ksp{u})) is the scattering potential defined in Ref. [118], and from which t{ui) can be calculated according 
to Eq. (1426). Notice that to perform such a calculation of t{uj) is not at all straight forward. The reason being 
that the reduced Rayleigh equation, through the Green's function Go(g, w), has poles at q = ±ksp{co). Details on the 
numerical method used for such calculations can be found in Ref. [118] and will therefore not be given here. 

As can be see from Fig. 21, ^ln|l + t(a;)|^^ is consistent, within the error bars, with a linear dependence on L. 

The solid line in this figure is a x^-fit to the numerical data. The slope of this curve is according to Eq. (144) related 
to the Anderson localization length, £{lo), as l/i{aj). Numerically the Anderson localization length for our system is 
found to be 

i{u)) = (5319 ±905) A. (147) 

This length should be compared to the one due to ohmic losses, which for our set of parameters is £e{oj) = 30.1 A. 

We have thus shown theoretically that by using specially designed surfaces, there might be hopes of observing the 
localization of surface plasmon polaritons at a randomly rough metal surface. 



C. Angular Intensity Correlations for the Scattered Light from Randomly Rough Surfaces 

It has been known for quite some time that when electromagnetic waves, all of the same frequency, are scattered 
from a random system, speckle patterns might be observed [78]. Such patterns are results of interference between 
waves scattered from diff'erent locations in the random medium. Prom studies off volume disordered systems, such 
patterns are known to contain a rather rich structure [78, 79]. In particular, it was predicted theoretically [79] for 
such scattering systems that there should exist three types of correlations — short-range correlations, long-range 
correlations, and infinitely-range correlations. These correlations were termed the C^^\ C^^', and C^^^-correlations, 
respectively, and they have all been observed experimentally [80-82] . 

In this section we will discuss speckle correlations, not for light scattered from volume disordered systems, but 
instead for light scattered from randomly rough surfaces. Examples of such speckle patterns obtained when an 
electromagnetic wave is scattered from a randomly rough surface are shown in Fig. 22. 

Let us start by considering a planar surface separating two different materials. Since the surface is planar, the 
scattering is completely understood as expressed through the celebrated Fresncl's formulae [10, 11]. Imagine an 
experiment where light is incident at an angle onto the interface. Since the surface is planar, all the light is scattered 
into the specular direction Og = Oq, and its intensity is given by Presnel formula. If we in a second experiment incident 
the light at an angle 0q = —Oq, i.e. at an angle that was the specular direction in the first experiment, all the light 
will be scattered into 0'^ = 9'q = —6q and its intensity is again given by Fresnel's formula. The scattered intensities 
in these two experiments are in fact equal. This is easily realized from the Fresnel formulae (866) by noting that the 
a-factors that they contain are unaffected by a change of sign in the momentum variables. Thus, if we know the result 
of the first experiment, say, we also know the outcome of the second one. In other words, these two intensities are 
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FIG. 22: Speckle patterns that result from the scattering of light of wavelength A = 612. 7nm incident on a rough silver surface 
at angles (a) 6*0 = 0° and (b) On — 10°. The Gaussian height-distributed surface was characterized by an RMS- height 5 = lOnm 
and a (Gaussian) correlation length a = 0.4yum. The length of the surface was L = 100/im and the dielectric constant of silver 
at the wavelength of the incident light is e(a;) = —17.2 + iO.498. 



perfectly correlated. We now introduce the momentum variables q and k related in the usual way to the scattering 
and incident angles respectively (see e.g. Eq. (149) below). Let the notation {q,k) denotes a corresponding pair 
of momenta variables where q is the scattered momentum and k the incident one. For our planar surface geometry 
we will thus have perfect correlation between the two scattering processes {q,k) and {q',k') if {q,k) = {—k',—q'). 
Furthermore, since any process, of course, is correlated with itself, we in addition will expect perfect correlation when 
{q,k) = {q' ,k'). 

The above example is rather trivial and well-known example of correlations in the scattered intensity from a planar 
surface. However, what happens to the intensity correlations if the surface is not planar, but instead randomly 
rough? This is an interesting and non-trivial question and we will address it in this section. In the discussion to 
be presented below we will be focusing on the angular correlations in the light scattered incoherently from a rough 
surface. Furthermore, we will mainly discuss the case where the surface is weakly rough. In particular we will try to 
answer the following question: When and under which conditions will the intensity scattered (incoherent) into the far 
field for different incident and scattering angles be related to one and other? 
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FIG. 23: The scattering system considered in the study of the angular correlation functions. 



1. Definition of the Angular Intensity Correlation Functions 

Let us start by introducing the unnormalized angular correlation function C(g, k\q' , k'), which we will define as^^ 

C{qMq'X) = {I{q\k)I{q'\k')) - {I{q\k)) {I{q'\k')) , (148) 

where I{q\k) denotes the intensity of the light scattered from the surface, and the angle brackets denote an average 
taken over an ensemble of realizations of the surface profile function C,{xi). Furthermore, the (lateral) momentum 
variables, q and fc, are both, in the radiative region (|g| < y/EoLo/c), understood to be related to the scattering and 
incident angles 9s and Oq respectively according to 



fc = VeoM — sin6'o, q = \J eo{uj) — sva.9s- (149) 

The primed momentum variables, q' and fc', arc related in a similar way to the primed angles 9'^ and 9'q. Theses 
angles, both primed and unprimed, are defined positive according to the convention indicated in Fig. 23. This figure 
also shows our scattering system consisting of a semi-infinite dielectric medium with a rough interface to vacuum. 
Furthermore, the intensity I{q\k) can be defined through the scattering matrix S{q\k) according to the formula 



im = ^n\sm\\ (150) 



Li V c 

where Li is the length of the xi-axis covered by the random surface. 

In many cases it is convenient to work with a normalized correlation function, E{q,k\q' ,k'), in contrast to the 
unnormalized one. The normalized angular intensity correlation function will we define by^^ 

"^^'^1*'''^ - {I{q\k)){I{q'\k')) • 

The lesson to be learned from the huge amount of research being conducted on correlation function in the field 
of random (bulk) disordered systems [79-81, 83, 84] is that there may exist correlations on many different length 
scales including short to infinite range correlations. Thus part of the challenge we are facing will be to separate these 
different contribution to C{q, k\q' , k') (or equivalently to S(g, k\q' , k')) from one another. 

The first step towards such a separation is to rewrite the correlation function in terms of the S-matrix. This is 
done by substituting the expression for the intensity, Eq. (150), into the definition of the correlation function and 
thus obtaining 

C{qMq',k') = [{\Siq\k)f |%'|fc')l') - {\S{q\kf) {\S{q'\k')f)] . (152) 



We have here suppressed any explicit reference to the polarization (the j^-index) since no confusion should result from doing so. All 
quantities in this section should be understood to be referring to one and the same polarization. 

It should be noticed that a somewhat different definition for the normalized angular intensity correlation function is used by some 
authors [89]. However, the advantage of the definition (151) is that it does not contain any (5-functions in the denominator. 
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Due to the stationarity of the surface profile function, the average of the S-matrix should be diagonal in q and k, 

{S{q\k)) = 2n5{q-k)S{k). (153) 
By now taking advantage of this relation in addition to the cumulant average [85, 86] 

{AB} = (AB) - (A) {B) , (154) 
the correlation function (152) can be written as 

C{q,k\q',k') = [\{SS{q\k)6S*{q'\k'))f + \{dS{q\k)6S{q'\k'))f 

+ {6S{q\k)SS*iq\k)6S{q'\k')6S*{q'\k')}] + s.t. (155a) 
where SS{q\k) denotes the incoherent component of the S-matrix defined as 

6S{q\k) = S{q\k) ~ {S{q\k)) . (155b) 

In Eq. (155a) the asterisks denote complex conjugate while s.t. means specular terms, i.e. terms that are proportional 
to S{q — k) and/or 5{q' — k'). Such terms will not be focused on here since we will concentrate on the incoherent part 
of the scattered light. With Eq. (155a) we can now write^^ 



C{q, k') = C^'\q, k\(/ , k') + C^^''\q, k\q' , k') + C^''\q, k\q' , k'), 



(156a) 



where 



and 



C^'\q,k\q',k') = ^'^\{6S{q\k)SS*{q'\k'))f, (156b) 
C^'''\q,k\q',k') = ^^'^ \{6S{q\k)5S{q'\k'))\\ (156c) 

2 

CW(9,fc|9',fc') = {SS{q\k)SS*{q\k)dS{q'\k')5S*{q'\k')}. (156d) 

Due to reasons which should be clear from the discussion below, the correlation functions in Eqs. (1566) and (156c) are 
termed short-range correlation functions, while the one in Eq. (156rf) contains contribution from long and infinite-range 
correlations. They will now be discussed in turn. 

2. Short Range Correlations for Weakly rough Surfaces 

In this subsection the short-range correlation functions, C^^^ and C^^"^ , for weakly rough surfaces will be discussed. 
These correlation functions are to leading order in the surface profile function a result of single scattering processes [91]. 
However, above leading order they will also receive contributions from multiple scattering. The long and infinite range 
correlations, C*^^-*, contain at least one multiple scattering process as wc will sec [91]. Therefore the "optical paths" 
involved in the processes leading to C^^' and C(i°) are typically shorter then those giving rise to C'^^K This is one 
of the reasons why the C^^^ and C^^^^ correlation functions are termed short-range correlation functions. Another 
reason stems from the fact that C*^"'^-' and C^^^^ arc both independent of the length of the random surface. In the next 
subsection we will demonstrate explicitly that the C'^'-correlation function is proportional to 1/Li. Hence, in the 
limit of a long surface the amplitude of the correlation function C'^' is neglectable compared to C'^^ and (7^^°^ 

The C^^^ Correlations Function; The Memory- and Reciprocal Memory-Effect 



Notice that equivalent expressions can be derived for the normalized correlation functions based on Eq. (151). 
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FIG. 24: Interpretation of the correlation condition for the short-range correlation functions C'^-* and C'^"'. The outgoing 
solid lines indicate the specular direction. The scattering process {q, k) that gives rise to the rnornenturn transfer Aqk = q — k 
might be correlated with the process {q', k') if Agk = Agi^' (C"^^') or if Agk = —Aqi^i (C^^""). 



At first sight, the expressions in Eqs. (156) might not seem too useful to us. However, they are as we now will 
try to explain. We will only be concerned about one-dimensional random surfaces, Q{xi), that are stationary and 
constitutes a Gaussian random process. Under this assumption the expression {5S{q\k)5S* {q'\k')) , contained in C^^\ 
will be proportional to a (5- function, i.e. 

{6S{q\k)5S*{q'\k')) oc 2-K5{q- k- q' + k'). (157) 

This is so due to the stationarity of the surface profile function C(2;i). To motivate this we recalling from Sect. Ill, 
or Ref. [55], that to lowest order in the surface profile function the scattering amplitude that is proportional to the 
S-matrix, is proportional to Q{q — k), where C denotes the Fourier transform of the surface profile function. Since 
{C{q)C* (k)) = 2TTS{q — k), Eq. (157), to lowest order, follows immediately. 

Thus, with Eq. (157) we find that the correlation fmiction can be written in the convenient form 

C^^)iq,k\q',k') = M(i^i^^:±^c«(g,fc|,',,'_, + fc). (158) 

Here Cq^^ is known as the envelope function of C^^^ and it is independent of the length Li of the surface. Notice that 
the C^^^ -correlation function can only be non- vanishing when the argument of the 5- function vanishes. Therefore, 
since 277(5(0) = Li, the (full) C*^^)-correlation function is also independent of the length of the surface. 

To see what the 5- function condition of Eq. (158) means physically, it is convenient to introduce the momentum 
transfer that can be associated with the scattering process. If the incident light has momentum k and the scattered 
light is described by the momentum variable q the momentum transfer is 

A,fc = q-k. (159) 

Such a scattering event we recall was earlier denoted by {q, k). Thus, what Eq. (158) says is that the two scattering 
processes {q, k) and (g', k') might have non-vanishing C^^^ correlations if and only if the two scattering events have 
the same momentum transfer, i.e. if and only if 

A,fc = Aga'. (160) 

This condition is depicted in Fig. 24. From the condition (160) it follows that if the incident momentum is changed 
from say k to k' = k + Ak, the entire speckle pattern shifts in such a way that any feature initial at q moves to 
q' = q + Aq. In terms of the angles of incidence and scattering, we have that if Oq is changed to O'q = 6o + AOq, any 
feature in the speckle pattern originally at Og is shifted to 6'^ = ds + Ads, where Ads = A9q{cos6q/ cos 6s) to first 
order in AOq. This effect can indeed be seen from the speckle patters presented in Figs. 22. 

It should in particular be noticed that condition (160) is satisfied if (i) k = k' and q = q' as well as if (ii) k = —q' 
and q = — fc'. These choices arc the ones mentioned in the beginning of this section for the scattering from a planar 
surface. It is interesting to notice that these for a planar surface trivial correlations, also holds true for the scattering 
from randomly rough surfaces, even though as should be noticed, their physical origin is rather different. Case (z) 
is kind of trivial since any scattering process should be perfectly correlated with itself. This effect is known in the 
literature as the memory-effect. Situation (ii), that doesn't seem that obvious at first, is, in fact, a consequence of 
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FIG. 25: The envelopes of the short-range correlation functions C^^'' (q, k\q' , k') (solid line) and C'^"^ (g, fc|g', fc') (dashed line) 
as a function of the scattering angle 9'^ for 9s = —10° and 6q — 20°. The angle 9'q is doterniincd from the 5- function constraint. 
The rough surface was a silver surface characterized by Gaussian height statistics of RMS-height S = 5nm. The correlation 
function was also Gaussian with a correlation length of a = lOOnm. The wavelength of the incident light was A = 457.9nm. At 
this wavelength the dielectric constant of silver is e(a;) = —7.5 -|- iO.24. (After Ref. [91].) 



the reciprocity of the S-matrix; S{q\k) = «S'(— fc| — q). Hence, when k = —q' and q = —k' there should be perfect 

correlations, and the effect is known as the reciprocal memory-effect. If the scattering system docs not possess any 
damping, the system also respects time-reversal symmetry. Due to this reason this latter effect is also known by some 
authors as the time-reversed memory effect. 

In Fig. 25 (solid line) we present the result of perturbative calculations [91] for the envelope of the C^^^ correlation 
function as a function of the scattering angle O',. for 9s = —10° and = 20°. The angle 9'q is determined from the 
5-function condition of Eq. (158). The incident wave was p-polarized, and the surface parameters are defined in the 
caption of this figure. Two well-pronounced peaks at scattering angles 9s = —20° and 9s = —10° are easily spotted 
in the envelope of C^^\ They corresponds respectively to the memory and reciprocal memory effect. It can in fact 
be shown that by instead considering the envelope of the normalized correlation function, one will have perfect 
correlation at the maximum point of these two peaks (see ee.g. Ref. [119]). 

Before continuing, we would like to point out that the memory and reciprocal memory effect seen in Fig. 25 are 
due to multiple scattering processes that involves surface plasmon polaritons. Thus, for an s-polarized wave incident 
onto a weakly rough metal surface, such peaks are not expected to be seen since in this case the incident wave cannot 
excite surface plasmon polariton at the rough surface [119]. However, for scattering of an ,s-polarized wave at a 
dielectric-dielectric interface the C^^^ may show peaks [120] even though no surface plasmon polaritons are involved. 
These peaks originate from multiple scattering processes involving so-called lateral waves [88] . 

Recently both the memory and reciprocal memory effect have been observed experimentally by West and 
O'Donnell [89] in the scattering of p-polarized light from a weakly rough, one-dimensional, random gold surface. 
We have reproduced one of their graphs in Fig. 26. 



The C^^°) Correlation functions 



We now focus on the C^^''-'-correlation function. This correlation function was originally overlooked in the early 
studies of correlation functions [87] due to the use of the factorization method [83]. By essentially duplicating the 
arguments used in arriving at Eq. (157), we find in an analogous way that 



{5S{q\k)5S{(f\k')) cx 2TT6{q-k + (f -k'), 



(161) 
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FIG. 26: Experimental measurements of the normalized correlation functions p\j{k,k' , /S.qk) (solid lines) and p^j{k, k' , Aqk) 
(dashed lines) as defined by West and O'DonncU [89] as function of the angle of incidence 9'^. Those correlation functions are 

these authors equivalent to our envelope functions S,l^\q,k\q' ,k') and sl^'^\q,k\q' ,k') (see Ref. [89] for details). The incident 
light had wavelength A = 6.12nm and the momentum transfer was Agk = 0.04(a;/c). The Gaussian height-distributed gold 
surface had RMS-height 5 ~ 15.5nm. Its correlation was characterized by a West-O'Donnell (rectangular) power spectrum of 
parameters k- = 0.83(ti;/c) and k+ = 1.30(a;/c). These values satisfy k- < kgp < k+ where ksp = 1.06(a;/c) is the surface 
plasmon polariton wave vector, and hence an incident wave should couple strongly to such modes. The memory and the 
time-revised memory peaks are indicated in these figures by A and TZ respectively. At these two positions we see that there 
are perfect correlations. (After Ref. [89].) 



with the consequence that we might write 

C^''\q,kW ,k') oc M(i^A±^:^C7r(9,fck',9' + 9 + fc). (162) 

Here C'^^^ {q, k\q', q' + q + k) is an envelope function, and both C^^"^ and its envelope Cq^°^ are independent of the 
length of the randomly rough surface. 

The presence of the ^-function on the right hand side of Eq. (162) is in terms of the momentum transfer equivalent 

to 

Agk = -Ag>k'. (163) 

What this condition implies for the speckle pattern is that if wc change the angle of incidence in such a way that k 
goes into k' = k + Ak, a feature originally at q = k — Aq will be shifted to q' = k' + Aq, i.e. to a point as much to one 
side of the new specular direction as the original point was on the other side of the original specular direction. For 
one and the same incident beam the C*'-^"^ correlation function therefore reflects the symmetry of the speckle pattern 
with respect to the specular direction (see Fig. 24). 
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The dashed Une in Fig. 25 shows the angular dependence, obtained from perturbation theory [91], for the envelope 
of C'^°). The parameters used to obtain these results were the same used to obtain the C^^^ correlation shown by the 
solid line in the same figure. It is seen that the Cq^^^ envelope is a smooth function of 6'^, and in particular does not 
show any peaks. Moreover, its amplitude is roughly of the same order of magnitude as the C^^"^ correlation function. 
This behavior is the same as the one found by West and O'Donnell [89] in their experimental investigation of the 
C^^°' envelope (Fig. 26). 

It should be pointed out that the C^^^^ correlation function has no known analogy within scattering from volume 
disordered system. This new type of correlations in surface scattering was first predicted from perturbation theory 
by Malyshkin et al. in 1997 [90, 91]. 

3. Long- and Infinite-Range Correlations 

We will now consider the last term of the right hand side of Eq. (156) that gives rise to C'^' . Due to the stationarity 
of the surface 

{SS{q\k)dS*{q\k)SS{q'\k')6S*{q'\k')} cx 2tt5{0) = Li. 
Hence, the correlation function itself, in light of Eq. (156(i), should behave as 

CW(g,fc|g',fc') oc ^. (164) 

It should be noticed that the C'^^^-correlation function is not constrained in its momentum variables through S- 
functions as wc saw earlier was the case for the short-range correlation functions. 

Even though we will not address this point explicitly here it has recently been shown that C'^^^ can be written as 
a sum of the three following terms [90, 91] 

C(^) {q, k\q',k') = C^l--^) {q, k\q\ k!) + C^^) (9, k\q' , k') + C^^^ (q, fcjq', k'). 

Here C'-^'^-' denotes a correlation function of intermediate-range, is a correlation function of long-range, while C*^^-* 
is an infinite-range correlation function. For explicit expressions for these three correlation functions the interested 
reader is directed to Rcfs. [90, 91] and [92]. In scattering from bulk disordered systems C^^^ [79, 81, 84] and C^^^ [79, 
80] have their analogies. However, the intermediate range correlation function, C^^'^\ predicted theoretically by 
Malyshkin et al. in 1997 [90, 91], has no equivalent in scattering from volume disordered systems. It is unique to 
scattering from randomly rough surfaces that support surface plasmon polaritons at the frequency of the incident 
light. An explicit example of such a scattering system is provided by a randomly rough metal surface in p-polarization. 

Based on a diagrammatic perturbation theoretical study, Malyshkin et al. [91] found that C^^'^^ shows a rather rich 
peak structure. Peaks in C^^'^^ are expected to occur for a number of cases in which a linear combination of three of 
the momenta q, k, q' and k' add up to ±ksp, where ksp, is the wave vector of the surface plasmon polariton. These 
condition arc summarized in Table II. In an expansion of C(i-5) in powers of the surface profile function the leading 
order is 0{(^). The intermediate-range correlation function C^^ "''^ is therefore for a weakly rough surface a result of 
correlations between a single scattering and a multiple scattering process that involves surface plasmon polaritons. 
In Fig. 27a we have plotted the intermediate-rang correlation function C^^'^^ for the randomly rough silver surface 
that lead to the results shown earlier in Fig. 25. In this graphs several peaks are easily seen. Their positions should 
be compared to the predictions that can be obtained from Table II. 

So far there is no experimental measurements for any of the correlations contained in C^^\ In fact such an 
experimental confirmation represents a real challenge to the experimentalists. The reason being that for long surface 
these correlations are rather small (due to Eq. (164)). In order to be able to observe them, one probably has to use a 
well-focused incident beam, or a short surface. 

Malyshkin et al. [91] also showed perturbatively that the C^^^ -correlation function should have a peak structure, 
while the infinite range correlation function, C'^^\ should be a smooth function of its arguments. This is seen from 
the perturbative results plotted in Fig. 27b (C^^^) and Fig. 27c (C^^^). The correlations described by the C^^^- 
correlations function are a result of correlation between two multiple scattering processes. For weakly rough metal 
surfaces this correlation function is dominated by double scattering processes. Its peaks arc associated with surface 
plasmon polaritons, as was found to be the case also for C^^'^-'. The peak conditions for C^^' are that two of the four 
momenta involved should add/subtract to zero. That is to say that for fixed k, q, and k', peaks are expected when 
q' = ±fc', q' = ztk or q' = ztq (see Table II). Also the infinite range correlations are due to multiple scattering events. 
What distinguish the long-range correlation, C*-^-* , from the infinite-range, C'^^^ , is that the latter involves at least one 
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FIG. 27: Perturbative results for the angular dependence of the correlation functions (a) C^^'^\q,k\q' ,k'), (b) C^'^\q,k\q' ,k') 
and (c) C^^\q, k\q' ,k') on the scattering angle for O3 = —10°, Oo = 20° and 6'q = 30°. The remaining parameters were the 
same used in Fig. 25. (After Ref. [91].) 



triple scattering process^^. For more details information about C^^'^\ C^'^\ and C^^\ the reader is invited to consult 
Refs. [90, 91] and [92]. 

4.. Angular Intensity Correlation Functions for Strongly Rough Surfaces 

Before closing this section, we would like to make a few remarks regarding strongly rough surfaces. Above we 
always assumed that the surface was a weakly rough metal surface. We saw that many of the interesting features of 
C{q. k\q' ,k') appeared due to excitations of surface plasmon polaritons. For strongly rough surfaces the excitation of 
surface plasmon polaritons, if any, is weak, and the dominating mechanism for multiple scattering from such surfaces 



The leading contribution to C(^) is of order in the surface profile function C{xi) [91]. 
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TABLE II: The peak conditions for the intermediate range C('-^) and long range correlation function C'^^ for a metallic 
one-dimensional surface. See text for details. 




-90.0 -60.0 -30.0 0.0 30.0 60.0 90.0 -90.0 -60.0 -30.0 0.0 30.0 60.0 90.0 



e's [deg] 9', [deg] 

FIG. 28: Rigorous numerical simulation results for the (a) Cg^^ and (b) Cq^""* envelopes as functions of 9's for 9o = 30° 
and 9s =0°. The angle 9'o was determined from the 5- function constraint. The s-polarized incident light had wavelength 
A = 632. 8nm. The randomly rough silver surface characterized by a (Gaussian) correlation length a = 3.85/Ltm. The RMS-height 

of the Gaussian height-distributed surface was 5 — 1.278/im (solid line) 5 — 0.1278/im (dashed line). As the RMS-height is 
increased one observes that the memory and reciprocal memory peaks start appearing in the envelope of C^^\ (After Ref. [119].) 



is multiple scattering of volume waves. As might have been guessed, multiple scattering of volume waves take over 
for strongly rough surfaces the role that surface plasmon polaritons had for weakly rough surfaces. These multiple 
scattered volume waves give rise to the memory and reciprocal memory effect for strongly rough surface. This is in 
fact the case for both p- and s-polarized incident light in contrast to what is the case for weakly rough surfaces. This is 
illustrated by the rigorous computer simulation results of Fig. 28a showing the Cq^^ envelope for s-polarized incident 
light [119]. It is seen from this figure that as the RMS-height is increased from a value corresponding to a weakly 

rough surface the memory and reciprocal memory peaks start to emerge in the Cq^^ envelope due to the increased 
contribution from multiple scattered volume waves. 

It should also be noticed, as was realized recently [119], that a measurement of the angular intensity correlations can 
provide valuable information regarding the statistical properties of the amplitude of the scattered field. In particular, it 
was shown that the short-range correlation function C^-^^^ is in a sense a measure of the non-circularity of the complex 
Gaussian statistics of the scattering matrix. If the random surface is such that only the C^^^ and C^^"^ correlation 
functions are observed, then S{q\k) obeys complex Gaussian statistics. If the random surface is such that only C^^^ 
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is observed, then S{q\k) obeys circular complex Gaussian statistics^^. This can indeed be seen from Fig. 28b, which 
shows that as the surface is made rougher, and therefore SS{q\k) approaches a circular complex Gaussian process, 
the C^^°'-correlation vanishes as compared to C^^-'. Finally, if the random surface is such that C^^'^\ C^^^ and C^^^ 
are observed in addition to both C^^^ and (7^^°^ then 5S{q\k) is not a Gaussian random process at all. However, 
which kind of statistics 6S{q\k) satisfies in this case is not clear for the moment. These results fits the findings from 
standard speckle theory [78, 93, 94] which assumes that the disorder is strong and that 6S{q\k) constitutes a circular 
complex Gaussian process. 

D. Second Harmonic Generation of Scattered Light 

So far in this section, we have discussed exclusively rough surface scattering phenomena that find their explanation 
within linear electromagnetic theory. There are still many exciting nonlinear [35] surface scattering effects that have 
to be addressed in the future. Such nonlinear studies are still at their early beginning. The studies that have been 
conducted so far on nonlinear surface scattering effects have mainly been related to the angular distribution of the 
scattered second harmonic generated light [10, 11]. In particular what have been studied are some new features in 
the backscattering directions of the second harmonic light. In this section we will discuss some of these results. The 
presentation given below follows closely the one given in Ref. [96] . 

It is well-known from solid state physics that an (infinite) homogeneous and isotropic metal has inversion symme- 
try [38, 116]. A consequence of this is that there is no nonlinear polarization in the bulk. If, however, the metal 
is semi-infinite with an interface to vacuum, say, the inversion symmetry is broken. Thus, a nonlinear polarization, 
different from zero, will exist close to the surface. As we move into the bulk of the metal, this effect will become 
smaller and smaller and finally vanish. Therefore, one might talk about a nonlinear surface layer which through 
nonlinear interactions will give rise to light that is scattered away from the rough surface at the second harmonic 
frequency. 

The scattering system that we will be considering is the by now standard one depicted in Fig. 6. This geometry is 
illuminated from the vacmim side, X3 > by a p-polarized planar wave of (fundamental) frequency lu. Only the 

p-polarized component of the scattered second harmonic generated light will be considered here, even though there 
also will exit a weak s-polarized component due to the nonlinear interaction at the surface. However, the p-polarized 
component represents the main contribution to the scattered light at the second harmonic frequency 2lo, and will 
therefore be our main concern here. Moreover it will be assumed that the generation of the second harmonic light 
does not influence the fleld at the fundamental frequency in any significant way. 

To motivate the study, we in Figs. 29 show some experimental results (open circles) due to K. A. O'Donnell and R. 
Torre [95] for the so-called normalized^^ intensity of the second harmonic light scattered incoherently from a strongly 
rough silver surface. The surface was characterized by Gaussian height statistics of RMS-height S = 1.81/xm and a 
Gaussian correlation function. The transverse correlation length was a = 3.4/im. The wavelength of the incident light 
was A = 27rc/a; = 1.064/um, while the angles of incidence considered were ^0 = 0°, ^0 = 6°, and ^0 = 15° as indicated 
in Fig. 29. 

The most noticeable feature of the experimental results (open circles) shown in Figs. 29 are, without question, the 
dips seen in the backscattering direction. It should be recalled that for the linear problem one gets at this scattering 
angle an enhanced backscattering peak (result not shown) similar to the one shown e.g in Fig. 13. So why do we have 
a dip for the second harmonic light at the backscattering direction, and not a peak? 

1. Strongly Rough Surfaces: A Numerical Simulation Approach to the Second Harmonic Generated Light 

Below we will with the help of numerical simulations try to get a deeper understanding of what causes these dips. 
The nonlinear layer existing along the surface is of microscopic dimensions. Since we are working with the macroscopic 
Maxwell's equations it is natural to assume that this layer is infinitely thin. Under this assumption, the effect of the 
nonlinear boundary layer is accounted for in the boundary conditions to be satisfied by the field, and its normal 



Two complex random variables A = Ai + iA2 and B = Bi + iB2 are said to be circular complex Gaussian if [93, 94] (AiBi) = {A2B2) 
and {A1B2) = - (A2B1). 

It can be shown that the total power scattered from a randomly rough surface at the second harmonic frequency is proportional to 
the square of the irrandiance, the incident power per unit area, on the surface. One therefore defines the normalized intensity of the 
scattered second harmonic light so that it is independent of the incident power. The analytic expressions for this quantity can be found 
in Eq. (34) of Ref. [96]. 
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FIG. 29: The mean normalized second harmonic intensity as a function of the scattering angle 6s for the scattering of p-polarizcd 
light from a randomly rough silver surface. The surface was characterized by a Gaussian height distribution of RMS-height 
5 = l.Sl/um, as well as a Gaussian correlation function of correlation length a = 3.4/xm. The dielectric constants were at 
the fundamental and second harmonic frequency e{u!) = —56.25 + i0.60 and e(2io) = —11.56 + 0.37 respectively. The thick 
lines represent the results of numerical simulations and the open circles represent the experimental results of O'DonncU and 
Torre [95]. The incident plane wave had a wavelength A — 1.064/im. In the numerical simulations the surface had length 
L = 40A and it was sampled with an interval Axi = A/20. The numerical results were averaging over = 2000 realizations 
of the surface, and the angles of incidence were (a) 6o = 0°, (b) 6o =6°, and (c) 9o = 15°. (After Ref. [96].). 



derivative, at the second harmonic frequency. These boundary conditions have jumps at the nonhnear interface, and 
their degree of discontinuity depends on the nonhnear polarization, or equivalently, on the parameters that describes 
this polarization. The form of the (nonhnear) boundary conditions at the second harmonic frequency 2u) can be shown 

to be [96] 

J-+(a;i|2u;)-^-(xi|2u;) = ^(.xi), (165a) 
J\f+ixi\2uj)-J\f-{xi\2u;) = B{xi), (165b) 

where the sources and Af^ have been defined in Eqs. (124). As before, the superscripts ± denote the sources 
evaluated just above (+) and below (— ) the rough surface defined by = ({xi). The functions A{xi) and B{xi) 
are related to the nonhnear polarization P(.Ti,a;3) through the integral of this quantity over the nonlinear boimdary 
layer [96]. To fully specify the nonlinearity of the problem, the polarization P(a;i,a;3) has to be specified. For instance 
for a free electron model, that we will consider here for simplicity, it takes on the form [97-99] 



P(a;i,a;3) = 7V (E • E) + /3E (V • E) . 



(166a) 
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Here the constants 7 and /? are defined as 



e%(rj_(xi,a;3)) 

7 = ^z:r7A ' (166b) 



= 3z!r:2' (166c) 



where no is the electron number density, rj_(xi, 2:3) is a vector normal to the local surface at point (xi, 0:3), and e and 
m are the charge and mass of the electron respectively. The explicit expressions, in this model, for A{xi) and B{xi) 
can be found in Ref. [96]. 

Since the surfaces used in the experiments leading to the results shown in Figs. 29 arc strongly rough, perturbation 
theory does not apply, and one has in theoretical studies to resort to rigorous numerical calculations of the second 
harmonic scattered light. Such kind of simulations are conducted on the basis of the rigorous simulation approach 
presented in Sect. Till. The calculations are now, however, made out of two main steps: First, one calculates the 
(linear) source functions JF^(xi|w) and A/'i/(xi|w); the field and it's normal derivative evaluated on the surface at the 
fundamental frequency uj. This is done exactly as described in Sect. Ill I. From the knowledge of the linear sources 
functions at the fundamental frequency, the right-hand-side of the boundary conditions (165) can be calculated since 
they depend directly on these source functions as well as on the form of the nonlinear polarization P(xi, 2:3) [96]. In 
all numerical results to be presented later in this section the form for the nonlinear polarization given by Eq. (166) 
will be used. With the functions A{xi) and B{xi) available, the nonlinear sources, J^^{xi\2lo) and M^{xi\2lo), are 
readily calculated from an approach similar to the one described in detail in Sect. Ill I. The only main difference is 
that now the boundary conditions to be used when coupling the two integral equations are the nonlinear boundary 
conditions given in Eqs. (165). With the source functions both for the fundamental and second harmonic frequency 
available, all interesting quantities about the scattering process, both linear and nonlinear, are easily obtained. The 
full details of this approach can be found in Ref. [96] . 

Based on this numerical approach, we compare in Figs. 29 the numerical simulation results (solid lines) obtained by 
Leyva-Luccro et al. [96] to the experimental results obtained by O'Donnell and Torre (open circles) [95]. The dielectric 
constants used in the simulations were at the fundamental frequency e{(jj) = — 56.25 + i0.60 and e{2u)) = — 11.56 + i0.37 
at the second harmonic frequency. Indeed by comparing the experimental and theoretical results shown in Figs. 29, a 
nice correspondence is observed both qualitatively and quantitatively. Particular in light of the oversimplified model 
used in the simulations for the nonlinear interaction, the agreement is no less then remarkable. 

From the experimental and theoretical results shown in Figs. 29a, a clear dip is seen in the incoherent component 
of the mean normalized second harmonic intensity for the backscattcring direction Og =0°. For the linear scattering 
problem, however, there is an enhancement at the same scattering angle. So what is the reason for the dip in the 
second harmonic light? O'Donnell and Torre [95], who conducted the experiments leading to the experimental results 
shown in Figs. 29, suggested that these dips were due to coherent effects. In particular they suggested that the dips 
originated from destructive interference between waves scattered multiple times in the valleys of the strongly rough 
surface. Since the numerical simulation approach seems to catch the main physics of the second harmonic generated 
light, it might therefore serve as a useful tool for testing the correctness of the suggestion made by of O'Donnell and 
Torre [95]. 

This can be done by applying a single scattering approximation to the generation of the second harmonic light. As 
described above, the numerical approach leading to the theoretical results shown as solid lines in Figs. 29, consists 
mainly of a linear and nonlinear stage where each stage is basically solved by some variant of the approach given 
in Sect. IIII. By using a single scattering approach, like the KirchhofF approximation [8, 106], to both stages of 
the calculation, a single scattering approximation for the full problem is obtained. The single scattering processes 
included in such a calculation is illustrated in Figs. 30. Notice that also unphysical scattering processes like the one 
shown in Fig. 30b, are included in this approximation. 

In Fig. 31 we present the consequence for the angular dependence of the normalized intensity {Ip{9 s^luj)) of only 
including single scattering processes in the second harmonic generation. From this figure it is easily seen that the 
intensity of the second harmonic generated light calculated in a single scattering approximation does not give rise to 
a dip (or peak) for the backscattcring direction. In fact the overall angular dependence of (/p(0s|2w)) in the single 
scattering approximation is quite different from the one obtained by the rigorous approach described above. Similar 
result holds for the other two angles incidence considered in Figs. 29. Hence, one may conclude that the dips present 
in the backscattcring direction of the incoherent component for the mean second harmonic generated light is not due 
to single scattering. It therefore has to be a multiple scattering phenomenon. 

To look more closely into this, the authors of Ref. [96] used an iterative approach for the linear part of the scattering 
problem which enabled them to calculate the scattered fields according to the order of the scattering process. Such 
a (Neumann-Liouville) iterative approach has been developed and used earlier in the literature [21, 100, 101]. For 
the nonlinear part of the calculation the rigorous simulation approach was used and thus all higher order scattering 
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FIG. 30: Diagrams illustrating two of the single scattering processes that produce the second harmonic scattered light in a 
single scattering approach. The double line black arrows represent light of frequency w, while the thick gray arrows represent 
light of frequency 2oj. Notice that the process in Fig. 30b is unphysical. 



processes were here taken into account. Some of the processes accounted for by this procedure and which give rise to 
the second harmonic hght is depicted in Figs. 32. 

We notice that the processes depicted in Figs. 32a and b represent single scattering in the linear part and are thus 
taken properly into account by using the standard Kirchhoff approximation [106] (for the linear part). However, for 
the paths shown in Figs. 32c and d, one needs to consider a pure double scattering approximation in order to include 
these processes properly. In Figs. 33 the simulation results for the second harmonic light (/p(^s|2a;)) are shown for the 
case where a single scattering (Fig. 33a) and a pure double scattering (Fig. 33b) approximation is used for the linear 
part of the scattering process. In both cases dips in the backscattcring direction are observed. In order to obtain the 
solid curve of Fig. 33c both single and double scattering processes were taken into account for the linear part of the 
calculation. This result would therefore include any interference effect between paths like those show in Figs. 32a-e. 
The dashed line in Fig. 33c is just the sum of the curves shown in Figs. 33a and b. It does therefore not contain any 
interference effects between type I paths (Figs. 32a b) and type II paths (Figs. 32c-d). That the two curves shown 
in Fig. 33c are so close to each other tells us that the interference between type I and type II paths are rather small 
(if any). Furthermore, paths of the type illustrated in Fig. 32e do not seem to be important, and they do not have 
coherent partners. 




FIG. 31: The mean normalized second harmonic intensity (^p(^s|2w))j„^ as a function of the scattering angle ds calculated in 
a single scattering approximation. The remaining parameters of the simulation were the same as those used in obtaining the 
results shown in Figs. 29. The angle of incidence was 6o =0°. (After Ref. [96].) 
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FIG. 32: Diagrams illustrating some of the multiple scattering processes that produce the second harmonic scattered light. 
The double line black arrows represent light of frequency w, while the thick gray arrows represent light of frequency 2a;. (After 
Ref. [96].). 



The numerical results presented so far seem to indicate that the behavior in the backscattering direction is affected 
by interference between the paths of either type I or type II. In the backscattering direction, there is no phase 
difference due to optical path difference between the two type I paths say. Similar argument hold for the type II 
paths. Hence any phase difference between the two paths has to come from phase shifts during the reflection. In the 
linear multiple-scattering processes giving rise to enhanced backscattering the phase shifts due to reflection will be the 
same for the two processes because the local Fresnel coefficients are even functions of the angle of incidence. Hence, 
the two paths in the backscattering direction will for the full linear problem both have the same phase and hence 
interferer constructively giving rise to the celebrated enhanced backscattering pack. However, for multiple scattering 
processes involving second harmonic generated light the situation is quite different. The reason for this is that the 
local nonlinear Presnel coefficient is not an even, but an odd function of the angle of incidence [96] . Hence, the phase 
difference between the two type I paths, say, will not be zero any more in general since the phases for these two 
paths will add instead of subtract. If this phase shift is positive in Fig. 32a, say, then it will be negative for the 
path shown in Fig. 32b since the local incident angles in the two cases have different signs and the local nonlinear 
Fresnel reflection coefficient is an odd function of the incident angle. Hence in the nonlinear case the phase difference 
in the backscattering direction is different from zero for the paths that seem to interfere. From the numerical results 
shown in this section they in fact seem to be close do tt out of phase resulting in destructive interference, or a dip as 
compared to its background at the backscattering direction. 



2. Weakly Rough Surfaces 

So far in this section we have presented both experimental and numerical results for the second harmonic gen- 
erated light scattered from strongly rough surfaces. There has also been conducted experiments for weakly rough 
surfaces [102]. The results are quite similar to the experimental results presented in Fig. 29. In particular, also for 
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FIG. 33: Calculations of the mean normalized second harmonic intensity as a function of the scattering angle 9s for the 
scattering of p-polarized light from a random silver surface where the linear part of the problem was solved by iteration. The 
incident angle of the light was = 0° and the other parameters of the simulation were as in Fig. 29. The curves have, (a) 
the single scattering contributions in the linear scattering and all contributions at the harmonic frequency, (b) pure double 
scattering contributions in the linear scattering and all contributions at the harmonic frequency, and (c) the single and double 
scattering contributions in the linear scattering and all contributions at the harmonic frequency. In (c), the curve shown with 
the dashed line represents the sum of the curves shown in (a) and (b). (After Ref. [96].) 



these weakly rough surfaces the second harmonic generated hght scattered diffusely showed a dip in the backscattering 
direction. However, in theoretical studies [103-105] both dips as well as peaks in the backscattering direction have 
been predicted. If it is a peak or dip depends on the values used for the nonlinear phenomenological constants. Even 
though predicted theoretically, only dips have so far been seen in experiments. 

For weakly rough surface the scattering processes giving rise to these dips (or peaks) are believed to be different 
for weakly and strongly rough surfaces. This situation resembles quite a bit the origin of the enhanced backscattering 
peak for weakly and strongly rough surfaces. Indeed, for weakly rough surfaces the origin of the dip in the intensity of 
the diffusely scattered light at frequency 2oj is intimately related to the excitation of surface plasmon polaritons at this 
frequency [104, 105]. Thus such dips are not to be expected for the second harmonic light generated in s-polarization 
from weakly rough surfaces. 



V. DIRECTIONS FOR FUTURE RESEARCH 



We have in the introduction to this review tried to give some glimpses of the many multiple scattering phenomena 
that may take place when electromagnetic waves are scattered by a randomly rough surface separating two media of 
different dielectric properties. Even though much is understood today when it comes to the rough surface scattering 
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problem, there are still, after a century of research efforts, many questions that have not been addressed and answered 
properly. Below we will therefore try to sketch out some directions for further research. 

We have exclusively considered one-dimensional surfaces. Naturally occurring surfaces arc mostly two-dimensional. 
Thus, the advance most needed in the field are techniques, either numerical or analytical, that accurately and fast 
can handle electromagnetic wave scattering from two-dimensional surfaces of varying RMS-height. Two-dimensional 
weakly rough surfaces can be treated by perturbation theory [108, 109], but if the surface is not weakly rough this 
approach is not adequate any more. In principle a general solution of the scattering problem can be formulated on the 
basis of a vector version of the extinction theorem [8, 11]. However, the resulting system of linear equations that needs 
to be solved in order to calculate the source functions is so big, and therefore require so much computer memory, that 
it for the moment is not practical in general [110]. Thus, one has to come up with new and more efficient methods for 
solving this kind of problems, or, the less appealing approach, to wait for advances in computer technology to make 
the extinction theorem approach tractable from a computational point of view. 

The scientific community dealing with wave scattering from disordered systems, seems to be divided into two 
separate groups: (i) those that deal with surface disordered systems and (ii) those that concentrate on systems with 
volume disorder. In the future these two "groups" have to be unified to a much higher degree then what is the case 
today in order to deal with scattering systems consisting of bulk disorder materials bounded by a random surface. 
Strictly speaking there has already been published some works for such "dual" disordered systems [111-114], but still 
more work, and in a more general framework, need to be done for such problems. 

An area that needs to be addressed further in the future is the inverse scattering problem[115] in contrast to the 
forward scattering problem that is the one that has received the most attention in theoretical studies of wave scattering 
from randomly rough surfaces. In the inverse scattering problem one has information, e.g. from experiments, about 
the angular dependence of the scattered light and one is interested in trying to reconstruct the surface profile function 
or its statistical properties. This problem is quite difficult and huge research efforts have been spent on it in related 
fields like remote sensing and seismic in order to try to find its solution. So far a general solution to the problem has 
not been found. 
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APPENDIX A: MATRIX ELEMENTS 



In this Appendix, some calculational details are presented for the matrix elements appearing in the matrix equa- 
tions (132) used to determine the source functions needed in the rigorous numerical simulation approach given in 
Sect. mi. 

Prom this section, Eqs. (133), we recall that these matrix elements are defined as 



■^tn = / dx\A±{U\x\), 

.AC/2 

= / du A±{£,rn\U+u), (Ala) 

■/-AC/2 
/■a-FA5/2 

Btvn = / dx'-^ -B±(^„|a;i), 

JS.^-A£,/2 
.AC/2 

= / duB±{U\^n + u) (Alb) 

J-AC/2 

where we in the last transition have made a change of variable u = xi — and where the kernels, according to 
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Eqs. (125) and (130), are given by 

^ (A2a) 



A±{xi\x\;u)) = lim 'y{x[)dn'G±{r\r' ■,ui) 
T)— >o+ 47r 



xa = C{xi) +r] 
4 = C(^') 



B±{xi\x[;oj) = lim -5- G±(r|r';a;) 
rj— >o+ 47r 



^ > (A2b) 



4 = C(^') 

with r = {xi,xs) and a similar expression holds for r', and 

e„ = -^+ A^, n= 1,2,3,. ..,7V, (A3) 

with = L/N. In the above expressions G±(r|r';a;) denote the free-space Green's functions for the Helmholtz 
equation. In 2-dimensions, as we will be considering here, it can be written as [56] 

G±(r|r» = iTTjyW (s±^^|r-r'|), (A4) 

where Hq^\z) is the Hankel function of the first kind and zeroth-order [39, 56]. By substituting this expression for 
the Green's function into Eqs. {A2) for the kernels, one gets 

, , , / ^ r f ^\ LO^ H['\x±ixi\x',)) 

-7-0+ V 4y c2 X±(a;i|a;'J 

X [{x, - x[)C{x[) - (C(a;i) - C{x[) + v)] , (A5a) 

B±{x,\x[;u;) = ^hm (^-'-^ Hi,%±{x,\x[)) , (A5b) 

where we have defined 

X±{xi\x[) = V^)'^ ^{x^-x[r + {C{x,)-ax[)+vy- (A5c) 

Notice that since the Hankel functions are divergent for vanishing argument [39, 56], so are the kernels A±{xi\x'i;u)) 

and B± {xi\x'i; u)) . However, fortunately these singularities are integrable, so the matrix elements ^^'^ ^mn 

fact non-singular everywhere and in particular when = ■ We will now show this and obtain explicit expressions 

for these matrix elements. 

We start by considering the off-diagonal elements where the kernels are non-singular. In this case, one may 
approximate the integrals in Eqs. (^1) by for example the midpoint method [40] with the result that (m ^ n) 

Ai„ = A±{U\U,oj), (A6a) 
B±„ = A^B±{UU,u), (A6b) 

where the expressions for the kernels are understood to be taken in the form Eqs. {A5). 

So now what about the diagonal elements where the kernels are singular? In order to calculate these elements, we 
start by noting that X±{£,m\^m + u), needed in order to evaluate the matrix element, can be written as 



= V^-l{U)\u\ + ... (A7) 



where we have Taylor expanded C(?m+'") and where we recall from Eq. (1206) that 7(a;i) = y^l -|- (^''^{xi). Furthermore, 
by advantage of the following (small argument) asymptotic expansions for the Hankel functions [39] 

i?«(z) = ^ +l + 0(22ln^), (A8a) 

= + - ln-+7+- --+0(^^ln^), (ASb) 
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where 7 = 0.5772157 ... is the Euler constant. 

With these expressions, it is now rather straight forward to derive, to obtain the matrix elements by integrating 
the resulting expressions term-by-term. To demonstrate this we start with the matrix element. With Eqs. (^47) 
and {A8a) and passing to limit ry — » 0+ whenever no singularities results and one gets 



B 



± 

mm 



rAi/2 



r.AC/2 



f 

= 2 

Jo 

|ln [ye±-l{im)u^ -h 7j 



2 2 



22- 

TT 



In 



/£±- 



c 2e 
c 2e 



(A9) 



Here in the last transition we have Eq. [ASa) one more. 

Furthermore, for the leading term of the diagonal elements of A one gets in a similar way from Eqs. (A86) and {Al) 



A?/2 

,2 



du A±{(,m\£,m + U) 



-e±^ lim / 



2i 



+ ... 



= lim 



2r, 



du 



+ 



ml J — 



du 



r,^o+ 2tt J_m - 2C{U)u + 1 

47r 72 

lim 

27r ?7^o+ 

1 

2 



^ jhn [tan-i(-C'(^™)+7(Cm)«)]H--l +^^^^"^^'"^ 



'47r72(e™) 



4^7' (Cm) 



To sum up we have for the matrix elements 



and 



I 4 

In these equations A±{^m\^n) and B±{^rn\^n) are given by Eqs. {A2) 



m n 
m = n. 



(AlO) 



(All) 



(A12) 



APPENDIX B: THE x-FUNCTIONS USED IN SMALL AMPLITUDE PERTURBATION THEORY 

In this appendix some of the lengthy formulae found in small amplitude perturbation theory. Sect. IIIF, are given. 
In particular we here give the first few x-functions found in Eqs. (88). We will now in the next two subsection 
explicitly give these functions for p and s-polarization. All explicit reference to the frequency u) has been suppressed. 
We have also for completeness used eq for the dielectric constant of the upper medium. In the case of vacuum this 
constant is £0 = 1. 
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1. P-polarization 



The three first functions in the set {Xp"^} are [46]: 

- eoa{q) 

2ao(fc) 



^ ^ (eq - e)^ a(g)ao(pi) + gpi 2ao(fc) [goa(pi)a(fc) - epik] 
£ao(g) + eoa(g) eao(pi) + eoQ:(pi) £ao(fc) + eoa{k) 

xf\(l\pAP2\k) = {2ea''{q)al{k) + [a'' {q) + al{k)\ [eoa{q)a{k) - eqk] 

-2eoqka{q)a{k)} 



sao{k) + eoa(fc) 



eao{q) + eoa{q) sctoipi) + soaipi) 

X {ea{pi) [al{k) - pik] + eoa{k) [a'^{pi) - pik] } 



eao{k) + eoa{k) 



£0 - £ 



£a!o(pi) + £oQ;(p 



— [a(pi)ao(P2) +PiP2]| 



sceo{p2) + £oa{p2) £Q;o(fc) + £oQ!(A;) ' 



(Bla) 



^° ^ [£oQ;(p2)a(fc) - £piA:] 77^T~^~T7T' (■^■'■'^) 



2. S-polarization 

Here follows the corresponding expressions for s-polarization [46] : 

y"HM) - -i— £2—^ ^°°''=) (B2al 

£0 — £ £o ~ ^ \ 2Q!o(fc) 



+3[a(9)-ao(pi)]+2- 



ao(pi) + Q:(pi) ao(p2) + Q:(p2) J ao(fc) + a{k) 

(B2c) 

With these expressions wc close this appendix. 
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